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ABSTRACT 


Triple  and  fourth-order  turbulence  products  obtained  from  coincident  fine-spatial-resolution 
three-velocity-component  laser-Doppler-anemometer  measurements  are  presented  for  a  three- 
dimensional  pressure-driven  turbulent  boundary  layer  from  y+  >  3  to  the  outer  edge  in  local  wall- 
shear-stress-direction  co-ordinates.  The  wind  tunnel  flow  is  produced  by  a  Ree  =  5940,  Ure,  =  27.5 
m/s,  upstream  turbulent  boundary  layer  on  a  flat  surface  that  approaches  a  wing-body  junction, 
although  no  results  near  the  junction  vortex  are  presented.  For  most  of  the  presented  locations, 
the  mean  velocities  and  stresses  from  this  data  set  were  previously  discussed  by  Olgmen  and 
Simpson  (1995a).  New  data  at  additional  locations  and  for  a  mean  two-dimensional  flow  are 
presented. 

The  flow  is  progressively  skewed  in  a  counterclockwise  direction  while  decelerating  and 
then  is  skewed  in  a  clockwise  direction  while  accelerating.  At  the  upstream  locations  the  flow 

gradient  angles  in  the  outer  layer  persist  in  the  upstream  direction.  Large  Uw,  -vw  are  produced, 

which  change  sign  at  the  downstream  locations.  The  Reynolds  shearing  stresses  lag  the  mean 
velocity  gradients. 

At  all  velocity  profile  stations,  in  the  nearest  wall  region  (y+  <  15)  sweeping  motions  (u  > 

0)  produce  substantial  positive  IF ,  Uv2 ,  uw2 ,  hFv  values.  As  the  upstream  near-wall  mean  flow 
is  skewed  progressively  in  a  counterclockwise  direction  relative  to  the  local  freestream,  wallward 

motions  of  outer  layer  fluid  produce  positive  TFw,  TFw,  ~FJvw,  w3;  generally  v3  is  negative  nearest 
the  wall.  At  later  stations  where  the  near-wall  flow  is  being  skewed  clockwise  back  toward  the 
original  direction,  TFw,  TFw,  -uvw,  vF  are  negative;  TF  is  positive  where  the  flow  is  being  skewed 
clockwise.  Positive  vwz  indicates  the  importance  of  motions  away  from  the  wall  (v  >  0).  Only  at 

stations  where  the  skewed  flow  is  returning  to  the  original  direction  is  vw 2 <  0  and  wallward 
motions  dominant.  In  the  outer  layer,  motions  away  from  the  wall  (v  >  0)  generally  account  for  the 
triple-product  behavior.  Selected  fourth-order  products  are  presented  which  influence  the  turbulent 
diffusion  of  turbulent  kinetic  energy  equation  triple  products. 

Results  are  plotted  in  6  different  co-ordinate  systems  in  Appendices  A-F. 
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CHAPTER  I.  INTRODUCTION 


Recent  studies  of  three-dimensional  turbulent  boundary  layers  (3D-TBL)  have  dealt  with 
many  issues  related  to  the  physics  of  turbulent  flows.  Reynolds-averaged  turbulence  closure 
models  for  3D-TBLs  are  still  a  main  issue.  The  anisotropy  of  the  turbulence  structure,  the 

variation  of  the  -I lv  and  -vw  shear  stresses  under  the  action  of  the  major  rates  of  strain 

dU/dy  and  dW/dy,  and  the  structural  parameters  relating  Reynolds  stresses  are  issues  that  have 
been  examined  (Olgmen  and  Simpson,  1993,  1995a).  Simplified  2-D  calculation  methods  that 
have  been  extended  to  3D-TBLs  do  not  seem  to  work  satisfactorily,  especially  with  adverse- 
pressure  gradients  and  strongly  skewed  flows.  Simplified  equations  relating  the  shear  stresses 
to  mean  flow  quantities  and/or  relating  higher  order  moments  to  lower  order  moments  in  an  ad 
hoc  fashion  do  not  seem  to  reflect  the  flow  physics  in  3D-TBLs  (Simpson,  1995).  Not  only  are 
more  data  of  stress  and  higher  order  moments  needed,  but  also  new  relations  among  important 
Reynolds-averaged  flow  variables  should  be  developed  from  data.  Measurements  in  3D-TBLs, 
including  higher  order  moments,  are  necessary  for  modeling  the  correct  nature  of  turbulence. 
Many  existing  data  sets  include  only  the  mean  velocity  and  Reynolds  stresses  since  the 
measurement  of  triple  or  higher  order  products  was  difficult.  Recent  data  sets  include  some  higher 
order  moment  information.  However,  the  very  near-wall  characteristics  are  not  reported  for  high 
momentum  thickness  Reynolds  numbers  (Ree)  due  to  physical  limitations  of  the  instrumentation 
or  are  available  for  mainly  low  Reynolds  number  flows. 

Bradshaw  and  Pontikos  (1985)  made  measurements  on  an  ‘infinite’  swept  wing 
configuration  similar  to  Elsenaar  and  Boelsma  flow  (1974)  at  4  locations  (Ree=5000).  Mean 
velocities,  Reynolds  stresses  and  triple  products  were  measured  with  single  hot-wire  and  cross¬ 
wire  hot-wire  probes  above  y+=230  and  the  skin  friction  measurements  were  made  with  a  surface 

fence.  The  shear  stress  angle  (SSA=arctan( vwfffl) )  lagged  behind  the  flow  gradient  angle 


(FGA=arctan(. 


dWfdy 

dU/dy 


)).  Townsend’s  A,  structural  parameter  (the  ratio  of  the 


magnitude  of  the 


shear  stress  to  the  twice  the  turbulent  kinetic  energy  A,  ={(-uv  f  +  (-vw  f  l(u2+v7+wz)  decreased 
with  the  effect  of  3-D  from  a  value  of  0.17  to  0.08  at  the  last  station.  Their  turbulent  kinetic  energy 
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(TKE)  balance  showed  that  turbulent  diffusion  increased  in  magnitude  in  3-D  flow  and  is  as  high 
as  the  advection  terms,  but  is  smaller  than  the  production  and  dissipation  terms.  The  flatness  of 

v  fluctuations  increased  towards  the  layer  outer  edge  with  values  between  15  to  20.  The-I/v 
stress  progressively  decreased  downstream;  possibly  due  to  a  tilting  of  the  2D  vortices  in  the  y-z 
plane  that  reduced  the  effectiveness  of  turbulent  mixing  (Bradshaw  and  Pontikos,  1985). 

Driver  and  Hebbar  (1 987)  studied  a  3-D  flow  generated  by  a  spinning  cylinder  (Wsurface=Ue) 
under  a  zero  pressure  gradient.  Measurements  were  made  on  a  stationary  section  of  the  cylinder 
downstream  of  the  spinning  section  using  a  3  component  LDV  at  Ree=6000.  Measurements  of 
the  mean  flow,  Reynolds  stresses,  and  triple  products  were  made  above  y+=40.  The  skin  friction 
was  measured  using  a  fence  gage.  Computation  of  the  flow  based  on  the  Reynolds  stress  model 

of  Launder-Reece-Rodi  (L-R-R)  (1975)  showed  that  the  model  overpredicted  the-uv  and  -vw 
stresses  throughout  the  layers.  The  eddy  viscosity  was  anisotropic.  Severely  underestimated  triple 
products  were  obtained  using  some  relations  of  Daly-Harlow  (1970)  and  L-R-R  (1975)  between 
triple  products  and  lower  order  products.  However,  the  pressure  rate-of-strain  model  of  Naot  et 
al.  (1973)  and  Launder-Reece-Rodi  (1975)  seemed  to  estimate  the  pressure  rate-of-strain  data 
variation  reasonably  well  throughout  the  layers.  Driver  and  Hebbar  estimated  the  pressure  rate-of- 
strain  data  from  the  stress  transport  equation  balance. 

Driver  and  Johnston  (1 990)  studied  both  pressure-gradient  and  three-dimensionality  effects 
on  turbulence  using  a  spinning  or  stationary  cylinder  with  zero  or  varying  streamwise  pressure 
gradients  generated  by  the  diverging  tunnel  geometry  and  controlled  with  boundary  layer  suction 
at  the  walls.  The  Ree  of  the  flow  varied  from  4000  (2D  case)  to  6000  (3D  case  D.S1). 
Downstream  of  the  spinning  cylinder,  measurements  of  the  mean  velocity,  Reynolds  stress  and 
triple  products  were  obtained  on  the  stationary  cylinder  on  which  the  three-dimensionality  decayed 
under  an  adverse  pressure  gradient.  Measurements  as  close  as  y+=8  were  obtained  with  a  3- 
velocity-component  LDV;  and  skin  friction  was  measured  with  an  oil-film  interferometer.  The  TKE 
balance  showed  that  the  production  and  dissipation  terms  were  of  the  same  order  of  magnitude 
and  convection  and  diffusion  terms  were  nearly  zero  (case  D.S1,  U  re(=  Wsurtace),  indicating  an 
equilibrium  flow.  The  value  of  A,  dropped  to  =0.11  in  the  outer  region  (case  D.S1)  and  near-wall 
values  were  as  low  as  0.05.  The  SSA  lags  the  FGA.  Launder-Reece-Rodi  (1975)  and  Naot- 
Shavit-Wolfshtein  (1973)  pressure-strain  models  were  able  to  estimate  the  pressure-strain  terms 
closely  for  that  flow.  However,  the  shear  stress  direction  was  estimated  poorly. 

Littel  and  Eaton  (1991)  studied  a  3-D  TBL  on  a  rotating  disk.  The  mean  velocity  was 
measured  using  a  3-hole  probe  and  turbulence  and  higher  order  statistics  were  measured  with 


cross-wires  and  single-wire  hot-wire  anemometry  at  several  different  Reynolds  numbers  at 
different  radial  locations  with  different  angular  velocity  of  the  rotating  disk.  For  radius=0.421 
meters,  Ree  was  between  2645  to  6065  and,  measurements  were  above  y+=1 12.  The  skin  friction 
was  calculated  using  a  log-law  fit.  The  value  of  A:  decreased  monotonically  away  from  the  wall. 
The  ratio  of  the  spanwise  eddy  viscosity  to  the  streamwise  eddy  viscosity  was  between  =0.8  and 
1 .1 ,  indicating  approximately  an  isotropic  eddy  viscosity.  The  TKE  balance  showed  that  turbulent 
diffusion  was  negligibly  small. 

Schwarz  and  Bradshaw  (1993,  1994)  studied  the  cross-flow  development  and  decay 
generated  by  a  30°  bend  at  Ree  =  4100.  Data  include  the  mean  velocity  obtained  with  a  three- 
hole  pressure  probe  and  Reynolds  stresses  and  triple  products  obtained  for  y/8  >  0.1  (y+=150) 
with  cross-wire  hot-wire  probe.  The  shear  stress  angle  lagged  the  flow  gradient  angle  and  the  A, 
structural  parameter  was  reduced  with  3-D  effects.  The  TKE  energy  balance  showed  that 
turbulent  diffusion  term  is  smaller  than  the  dissipation  and  production  terms. 

Flack  and  Johnston  (1 993)  studied  the  3-D  TBLs  produced  by  flow  in  a  30°  bend  and  flow 
approaching  a  swept  45°  forward  facing  step  in  a  water  tunnel.  The  Ree  varied  from  =  1300  to 
1700  for  the  30°  bend  case  and  =  1400  to  1850  for  the  forward  facing  step  case.  Skin  friction  at 
the  wall  was  calculated  using  the  total  velocity  gradient  and  the  dynamic  viscosity  in  the  sublayer. 
Measurements  of  the  mean  velocity,  Reynolds  stresses  and  ten  triple  products  were  made  with 
a  3-velocity-component  LDV  as  close  as  y+=1.5.  The  SSA  lagged  the  FGA  in  the  outer  layer, 
while  below  y+=10  the  SSA  leads  the  FGA.  The  A,  parameter  decreased  with  increased  3D 
effects,  although  both  the  TKE  and  shear  stress  magnitude  increased.  The  distribution  of  some 
of  the  triple  products  in  the  layers  presented  in  that  study  are  qualitatively  similar  to  the  quantities 
presented  in  this  study. 

Here  higher-order  moment  results  are  presented  and  discussed  for  a  3-D  pressure-driven 
turbulent  boundary  layer  generated  by  a  3:2  elliptic-nose/NACA  0020  tail  wing  (chord  =  30  cm, 
thickness  =  7.17  cm,  height  =  22.5  cm)  perpendicular  to  a  flat  plate.  The  Reynolds  number  based 
on  momentum  thickness  measured  0.75  chord  upstream  of  the  wing  on  the  centerline  is  5940  for 
a  nominal  air  speed  of  27.5  m/sec.  LDV  data  were  obtained  at  9  stations  along  a  line  defined  by 

translating  along  the  flow-angle  direction  at  the  y  locations  in  each  profile  where  the  TF  normal 
stress  is  maximum.  This  location  was  conjectured  as  the  y  location  where  the  inner  and  outer 
layer  turbulence  structure  overlaps.  The  line  was  defined  with  hot-wire  measurements  at  18 
stations.  The  flow  studied  is  a  unique  data  set  in  which  the  lateral  pressure  gradients  that  turn 
the  flow  change  their  sign  and  thus  turn  the  very  near-wall  flow  back  toward  its  initial  direction. 
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Figure  1  shows  a  portion  of  the  wing,  the  measurement  locations,  and  the  local  free- 
stream  xFS  and  the  local  wall-stress  coordinate  xwc  directions  at  each  location.  The  xFS  coordinate 
is  in  the  direction  of  the  local  inviscid  velocity  immediately  outside  the  boundary  layer.  The  xwc 
coordinate  is  in  the  direction  of  the  flow  angle  at  the  first  nearest-wall  measurement  location  of 
each  station  profile.  The  y  coordinate  is  normal  to  the  surface  in  both  coordinates  and  the  zFS  and 
zwc  coordinates  form  right-handed  systems.  This  figure  shows  that  at  stations  6  to  9  the  wall 
coordinates  and  local  free-stream  coordinates  are  very  close  to  each  other  (table  1).  Mean 
velocities,  Reynolds  stresses,  and  skin  friction  values  as  well  as  the  flow  angle,  flow  gradient 
angle  and  shear  stress  angle  at  stations  1  through  7  were  previously  discussed  by  Olgmen  and 
Simpson  (1995a)  in  local  free-stream  coordinates.  For  completeness,  for  reference  and  for  some 
new  insights,  mean  velocity  components  (Section  II),  SSA,  FGA,  and  FA  (Section  III)  and  the 
Reynolds  stresses  (Section  IV)  are  presented  here  in  wall-stress  coordinates.  Besides  these 
previously  discussed  data,  2  new  additional  LDV  profiles  at  stations  8  and  9  are  also  presented. 
The  new  data  contain  the  mean  velocities,  Reynolds’  stresses,  10  triple  products  and  15 
quadruple  products,  where  the  3-D  flow  is  relaxing  downstream  of  station  7.  A  2-D  flow  profile 
was  also  measured  with  the  LDV  in  the  vicinity  of  the  station  5  without  the  wing  present  in  order 
to  document  the  variation  of  higher  order  terms  near  the  wall.  Data  presented  here  satisfy  the 
realizability  conditions  described  by  Schumann  (1977).  All  quantities  except  the  mean  velocities 
and  normal  stresses  are  nondimensionalized  with  respective  powers  of  the  normal  stresses  in 
correlation  coefficient  form.  Mean  velocity  components  and  normal  stresses  are 


nondimensionalized  by  respective  powers  of  skin-friction  velocity=u,  (ut= 


,tMr=wall  shear 


stress  magnitude,  p=density).  All  quantities  are  plotted  vs.  y+  (y*=—,  y=distance  from  the  wall 

v 

v=kinematic  viscosity).  At  stations  1  through  7  the  uT  skin  friction  velocity  was  obtained  with  an 
oil-film  interferometry  technique  by  Ailinger  (1990).  At  stations  8  and  9  skin  friction  is  calculated 
using  the  near  wall  mean  velocity  gradient  and  the  dynamic  viscosity.  For  the  2-D  profile  Coles 
law-of-the-wall  profile  with  constants  k=0.41  and  C=5  was  used  to  calculate  uT.  The  uncertainties 
in  the  measured  quantities  with  20:1  odds  are  given  in  table  2.  Measurements  were  made  with 
a  3-velocity-component  LDV  system  which  is  a  sub-system  of  a  5-velocity-component  LDV  system 
designed  to  study  near-wall  turbulence  (Olgmen  and  Simpson,  1995  b).  The  probe  volume 
diameter  was  30  pm  and  measurements  were  made  as  close  as  y+=3  to  the  wall  with  a  (+20, -5|i) 


4 


uncertainty  in  y  location.  The  surface  pressure  distribution  on  the  wall  was  measured  at  closely 
spaced  locations.  An  oil-flow  visualization  technique  was  used  to  reveal  the  skin-friction  line 
characteristics  for  this  flow.  Skin-friction  coefficient  magnitude  and  direction  were  reported  from 
measurements  by  an  oil-film  interferometry  technique.  However,  due  to  the  direction  uncertainty, 
the  wall-stress  direction  was  taken  as  the  mean  flow  direction  of  the  closest  to  wall  point  in  the 
profiles. 

For  reference,  the  conclusions  of  the  previous  work  (Olgmen  and  Simpson,  1995a)  are 
summarized  here.  The  flow  is  highly  anisotropic,  with  an  anisotropy  ratio  (ratio  of  the  lateral  to 
longitudinal  eddy  viscosities)  different  than  1.  In  wall-stress  direction  coordinates,  the  anisotropy 
constant  is  close  to  0.6,  although  it  varies  through  the  layers  and  from  station  to  station.  The 
shear  stress  angle  lags  behind  the  flow  gradient  angle.  Townsend’s  structural  parameter  A1  is 
observed  to  be  less  than  0.15  at  the  outer  regions  of  the  layers,  and  substantially  less  than  0.15 
near  the  wall,  with  double-peaked  profiles  at  some  stations.  The  ratio  of  the  vertical  normal  stress 

to  the  shear  stress  magnitude  (v^lx/p  | )  is  about  1 .6  in  the  y+=  1 00  to  1 000  range.  Also,  the  ratio 

of  the  shear  stresses  -vw/-Tiw  varies  semi-logarithmically  and  the  ratio  of  the  production  of 

these  stresses  C=  -vwproduction  lJuwproduction  is  about  constant  throughout  the  layers,  although  it 
varies  from  station  to  station.  For  further  details  see  Olgmen  and  Simpson  (1995a). 


CHAPTER  II.  MEAN  FLOW  AND  REYNOLDS  STRESS 
RESULTS 


The  U/uT  profiles  at  all  stations  follow  each  other  closely  up  to  y+=9  (figure  2).  Above  this 
height,  station  1,  2  and  3  values  closely  follow  each  other  up  to  y+=100.  These  are  the  stations 
where  three-dimensionality  of  the  flow  starts  to  develop.  The  U//ut  profiles  overlap  each  other  with 
a  semi-logarithmic  variation  in  the  y+=30-150  range,  suggesting  a  similarity  relation.  At 
downstream  stations  there  is  no  similarity.  Starting  with  station  4  the  profiles  become  fuller  and 
the  semi-logarithmic  regions  of  the  profiles  shift  to  higher  locations  in  the  profiles,  especially  at 
stations  8  and  9  this  effect  is  more  pronounced. 

This  presentation  of  the  data  in  wall  coordinates  is  in  the  law-of-the-wall  similarity  variables 
of  Coles  (1956)  and  van  den  Berg  (1975).  Nine  similarity  relations  including  these  two  relations 
were  examined  with  nine  different  data  sets  and  no  strong  evidence  for  a  law-of-the-wall  similarity 
relation  could  be  found  (Oilmen  and  Simpson,  1992).  Among  these  relations,  the  Johnston  (1960) 
similarity  relation  for  the  component  Q  =UFS/ cos(plv-pFS)  in  the  local  free-stream  direction  seemed 
to  fit  those  data  better  than  the  others: 


Q*  = 


U, 


FS 


- - - =B.  In(^)+S2 

ux  cos(p„-pra)  v 


(1) 


Johnston’s  relation  also  fits  the  present  data  best  (figure  3)  for  the  first  3  stations  with  B,  and  B2 
constants  close  to  the  Coles  and  Hirst  (1968)  constants  of  k  =  0.41  (=1/Bt)  and  C  =  5.0(=B2).  This 
suggests  that  a  law-of-the-wall  semi-logarithmic  profile  with  universal  constants  only  exists  for  the 
initial  stages  of  a  mean  3-D  flow. 

The  V/u,  profile  at  the  first  station  gradually  increases  from  the  wall  to  the  layer  edge 
(figure  4).  At  station  2  and  downstream  the  maxima  of  the  profiles  are  within  y+=600  to  1000  and 
grow  with  increasing  3-D  effects  as  the  local  free-stream  flow  accelerates.  Near  the  layer  edge 
at  each  station  the  values  decrease.  The  scale  of  the  V/uT  is  1/20th  of  U/uT  plots  and  in  the  region 
below  y+  <  100  the  mean  velocity  value  is  nearly  zero  for  stations  1  through  7.  This  shows  that 
below  y+=  100  the  vertical  motion  of  the  flow  is  due  to  turbulent  fluctuations.  At  stations  1  through 
8  up  to  y+=300  the  magnitude  of  V/uT  is  lower  than  the  2-D  value.  At  station  9  the  magnitude  is 
twice  of  the  magnitude  at  station  8  due  to  the  effect  of  the  horse-shoe  vortex  (Devenport  and 


Simpson,  1992)  lifting  the  flow  away  from  the  wall. 

Due  to  the  definition  of  the  wall-stress  coordinates,  W/ut  profiles  are  close  to  each  other 
very  near  the  wall  with  a  magnitude  of  zero  (figure  5).  At  the  first  3  stations  W/ut  monotonically 
increases  from  the  wall  to  the  layer  edge  since  the  flow  in  the  entire  layer  turns  in  one  direction. 
Between  stations  3  and  4  the  lateral  pressure  gradient  in  wall-stress  coordinates  changes  sign. 
At  station  4  the  W/ut  values  are  close  to  zero  up  to  y+=40.  Above  this  y  location,  values 
monotonically  increase.  At  stations  further  downstream  the  effect  of  the  sign  change  of  the  lateral 
pressure  gradient  is  felt  most  near  the  wall.  This  results  in  negative  W/ut  values.  The  magnitude 
of  the  absolute  maximum  progressively  increases  downstream  and  maximum  points  in  the  profiles 
shift  upward  in  the  layers.  The  Reynolds-averaged  Navier-Stokes  equations  show  that  the  mean 
flow  is  immediately  affected  by  the  mean  pressure  gradients.  The  pressure  force  is  most  effective 
on  the  near-wall  flow  where  the  momentum  of  the  flow  is  lowest.  The  lateral  pressure  gradient 
in  wall-coordinates  increases  in  magnitude  proceeding  downstream  of  station  4  till  station  7,  and 
thus  the  effect  of  the  pressure  gradient  is  felt  in  a  larger  y  range  near  the  wall.  At  stations  8  and 
9  the  lateral  pressure  gradient  is  negative  with  smaller  magnitudes.  Some  length  scales  calculated 
using  the  mean  velocities  are  given  in  table  3. 

II  A.  FLOW  ANGLE,  FLOW  GRADIENT  ANGLE,  SHEAR  STRESS 
ANGLE 

The  flow  angle,  flow  gradient  angle  and  the  shear  stress  angle  are  defined  as 

FA=arctan(J^)  ;  FGA  =arctan(  ;  SSA=arctan(J^)  (2) 

U  aU/ay  jjy 

and  are  plotted  against  y+  in  wall  stress  coordinates  (figure  6).  At  the  first  3  stations  the  flow 
angles  monotonically  increase  towards  the  layer  edge.  At  station  4  the  flow  angle  is  very  close 
to  zero  till  y+=40.  At  downstream  stations  the  flow  angles  are  negative  near  the  wall  due  to  the 
crossflow  pressure  gradient  effect  and  the  absolute  maximum  value  of  the  profiles  increases  till 
station  6.  Downstream  of  station  6  absolute  flow  angle  maxima  decrease.  The  location  of  the 
absolute  maximum  progressively  shifts  upward  in  the  layers.  The  flow  angle  at  the  layer  edge  in 
wall-stress  coordinates  increases  till  station  3.  At  downstream  stations  the  layer  edge  flow  angle 
decreases  till  station  9.  At  stations  6  through  9  the  wall-stress  coordinates  and  the  free-stream 
coordinates  are  oriented  in  close  proximity  to  one  another.  This  results  in  a  unique  flow  feature. 
While  the  near-wall  flow  and  the  free-stream  flow  are  in  directions  close  to  one  another,  the 


intralayer  flow  is  in  different  directions. 

At  the  first  2  stations  the  flow  gradient  angle  is  approximately  a  constant  above  y+=100. 
Further  downstream  they  show  a  maximum  peak  at  increasing  distances  away  from  the  wall.  At 
stations  5  to  9  minima  are  also  observed.  For  the  first  4  stations  the  shear  stress  angle  above 
y+=100  is  about  constant  in  the  tunnel  co-ordinate  direction  (Simpson,  1995).  At  downstream 
stations,  above  this  y+  location  the  shear  stress  angle  profiles  form  peaks.  At  every  station,  the 
shear-stress  angle  lags  behind  the  flow  gradient  angle.  Very  near  the  wall  large  differences  are 
observed  for  the  downstream  stations.  At  stations  8  and  9  the  difference  is  highly  increased.  At 
station  9  SSA  is  approximately  at  -90°  with  respect  to  the  FGA.  The  shear  stress  angle  responds 
to  the  changes  in  flow  gradient  angle  in  the  same  direction,  but  with  a  smaller  change  in 
magnitude.  At  station  5  through  7  this  results  in  the  intersection  of  the  shear  stress  angle  and 
flow  gradient  angle  profiles  at  2  points,  which  is  discussed  later  in  the  text.  At  station  8  the  SSA 
and  FGA  follow  each  other  closely  in  a  region  y+=25  to  300.  At  station  9  there  exists  a  similar 
region  above  y+=100  till  450,  where  the  main  stress  is  in  the  direction  of  the  main  strain. 

II  B.  REYNOLDS  STRESSES 

The  77V  normal  stress  profiles  steeply  increase  from  the  wall  to  their  peak  at  y+=15  to 
20  and  then  decrease  up  to  y+=60  to  450  at  different  stations  (figure  7).  The  fluctuation  continuity 
equation  dictates  that  nearest  the  wall  u  fluctuating  velocity  varies  as  y.  The  plateau  region  where 77V2 

is  approximately  constant  varies  at  different  stations.  The  maximum  in  the  profiles  decreases  till 
station  4  due  to  adverse  pressure  gradient  effects  and  then  increases  back  to  station  1  values 
downstream  due  to  the  accelerating  flow.  However,  above  y+=100  the  values  decrease 
monotonically  downstream  of  station  1  till  station  7.  The  decrease  at  station  7  is  approximately 

60%  of  the  station  1  value.  This  suggests  that  there  are  two  zones  in  the  flow  where  7T7  is 
distinctly  different:  1)  below  y+=100  which  is  highly  affected  by  the  pressure  gradients  and  wall 
proximity,  and  2)  above  y+=100  where  the  flow  relaxes  with  more  pronounced  3-D  effects.  The 
streamwise  pressure  gradient  is  favorable  downstream  of  station  2  until  past  station  6.  Between 

stations  6  and  7  ( dCp/dx)wc  changes  sign.  At  downstream  stations  77V2  increases  both  in  the 

inner  and  the  outer  regions.  The  peak  location  at  each  profile  is  between  y+=15  to  20.  The 
maximum  values  for  1 ,  6,  and  7  are  close  to  the  peak  value  for  2-D  flow.  At  stations  2  through 
5  the  peak  value  is  lower  and  at  stations  8  and  9  the  peak  value  is  higher  than  2-D  values.  Above 
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y+=200,  the  7F/u? values  are  lower  than  the  2-D  values  except  at  stations  1,  2  and  3. 


Profiles  of  v*/u?  monotonically  increase  to  a  maximum  between  y+=300  to  400  and  then 


decrease  toward  the  freestream  (figure  8).  Downstream  of  station  1  up  to  station  5  the  maximum 
value  decreases.  The  station  6  maximum  is  higher  than  the  maxima  at  stations  5  and  7.  At 
stations  8  and  9  the  profiles  show  two  maxima:  one  around  y+=60  and  one  around  y+=400.  The 
near  wall  values  at  these  stations  are  higher  than  station  1  values  below  y+=80,  and  away  from 


the  wall  the  values  are  in  between  the  station  3  and  4  values.  The  v 2/u?  profiles  of  stations  4 

through  7  are  lower  compared  to  the  2-D  values  and  other  stations  have  higher  values  compared 
to  2-D  values  throughout  the  layers. 


The  w^/uf  profiles  reach  maxima  at  a  lower  y+  location  than  the  v^/uf  stress  profiles 

(figure  9).  Values  at  station  2  are  higher  with  respect  to  station  1  and  station  3  values  are  close 
to  station  2  values.  Proceeding  downstream  up  to  station  7  values  decrease.  At  stations  8  and 
9  the  values  closely  follow  each  other  and  are  increased  close  to  station  4  values  above  y+=100, 
below  y*~30  the  values  are  close  to  station  2  values.  Compared  to  2-D  profile  all  stations  except 
station  7  have  values  higher  than  2-D  values  above  y+=1 00.  However  near  the  wall  at  stations 

4  through  7,  Wlu?\s  decreased  compared  to  the  2-D  profile,  presumably  due  to  the  increased 


favorable  pressure  gradient’s  effect  on  the  fluctuating  velocities. 

The  Reynolds-stress  transport  equation  with  summation-subscript  notation  can  be  written 
as  (Daly,  and  Harlow,  1970): 

dufi  sup,  — du.  — at/ 

dt  rlx^  '  lx,  1  ^ 
convection  production 


redistribution 


_8  (ufiu,) 
dx, 

turbulent  diffusion 


_1  MPUj)  J(pu,) 
p  dx,  dXj 
pressure  diffusion 


tfUju, 

+v _ — 

dx? 

viscous  diffusion 


„  du,  dU, 

-2v — 1 _ L 

dx,  dx, 

viscous  dissipation 


(3) 


The  production  rate  terms  for  the  IF,  v7  and'w2  stresses  can  be  written  as  (-2 uv^L), 

dy 

(~2v2^X.),  and  (-2 vw^.),  assuming  that  the  variations  in  the  x  and  z  directions  are  much 
dy  dy 

smaller  than  the  variations  in  y  direction.  The  relative  magnitudes  of  these  mean  velocity 

gradients  are  |^|>|A!l^|>|^|  and  the  stresses  are  v2> I Uv\> i vw\ .  This  results  in  th eF 
dy  dy  dy 

production  rate  being  much  higher  than  that  for  v2  and  w2.  The  equations  for  the  production 

terms  show  that  for  the  type  of  flow  studied  here  while  the  production  of  IF  is  never  zero  in  the 

layers  the  w2  production  is  zero  at  locations  where  vvv=0  and/or  dWIdy  is  0.  These  locations 

correspond  to  the  locations  where  the  shear  stress  angle  is  zero  ( arctan(vw/uv)  =0)  and/or  the 

flow  gradient  angle  is  zero  (arctan(..^^)=0).  Production  of  v2  is  zero  at  locations  where 

dU/dy 

(dV/dtf  =0.  At  stations  2  and  downstream  stations  this  location  is  below  y+=8.  The  production  oiv2 
is  an  order  of  magnitude  smaller  than  that  of  w2  and  two  orders  of  magnitude  smaller  than  that 
of  IF.  The  only  other  stress  transport  equation  term  reported  here  is  for  the  turbulent  diffusion. 

Since  d{UjJp,)ldx,  with  I  =  2  (in  the  perpendicular  direction)  is  the  main  turbulent  diffusion  term, 
where  it  is  zero  corresponds  to  the  locations  where  the  turbulent  diffusion  is  approximately  zero. 

At  all  stations  the  Uvl (IFF2) 1/2  correlation  coefficient  show  double  peaks;  one  near  the 
wall  below  y+=50  and  one  around  y+=300  to  1000  range  (figure  10).  The  near  wall  magnitude 
decreases  till  station  3.  Stations  4  and  5  show  the  magnitude  increase  due  to  the  favorable 
pressure  gradient.  Station  6  has  values  close  to  station  5.  Station  7  has  the  highest  near-wall 
magnitude.  At  stations  8  and  9  the  correlation  magnitude  decreases  again.  The 

Uvl  (IF  v2)'12  shows  the  combined  effect  of  streamwise  pressure  gradient  on  shear  stress  and 
normal  stresses.  The  2-D  data  values  reach  at  a  plateau  value  of  0.3  around  y+=30;  a  value  lower 
than  vast  amounts  of  previous  data  on  2-D  turbulent  boundary  layers.  To  resolve  this  ambiguity 
four  different  profiles  were  taken  which  did  not  change  the  results  presented.  Re-reduction  of  the 
raw  data  without  cleaning  the  noise  at  the  skirts  of  individual  histograms  of  U,  V,  W  velocity  data 
and  without  using  the  velocity  bias  weighting  scheme  resulted  in  lower  correlation  values.  Crude 


rotation  of  the  raw  data  obtained  in  optical  coordinates  to  tunnel  coordinates  without  using  the 
accurate  transformation  scheme  also  did  not  have  any  effect  on  the  results.  Also,  the  coincidence 
window  of  the  frequency  processors  were  varied  between  0.8  psec  to  10  psec  without  any 
changes  in  the  measured  quantities. 

At  the  2-D  station  I Iv/u?  is  approximately  a  constant  between  y+=  40  to  200  with  a  value 

of  -0.9  which  is  a  rather  high  value  than  -1 .0  obtained  in  the  literature.  At  station  1  uv/u?  reaches 
a  plateau  region  at  y+=40  at  about  -1 .2  and  stays  constant  up  to  y+=400.  An  effect  of  an  adverse 
pressure  gradient  on  2-D  TBLs  is  to  reduce  -uv  near  the  wall  and  increase  it  in  the  outer  region 

(Simpson  et  al.,  1981).  Above  y+=100,  -  uv/uf  increases  as  the  flow  moves  downstream.  However 

below  y  =*100,  -  uv/ i/T  increases  up  to  station  4,  and  is  followed  by  a  decrease  at  downstream 

stations.  This  decrease  below  y+=1 00  and  the  increase  above  this  y  location,  accompanied  with 
the  fact  that  the  stresses  both  at  the  wall  and  at  the  layer  edge  must  be  zero,  result  in  double- 
peaked  stress  distributions  at  stations  4  through  9.  The  minimum  value  obtained  is  at  station  8 
with  a  value  of  =-1.5  at  y+=35.  A  local  minimum  at  y+=250  with  values  close  to  that  of  station  7 
is  seen.  A  local  maximum  is  observed  at  y+=550  at  the  same  station.  At  station  9  the  peak 
minimum  value  is  decreased,  and  between  y+=300  to  1000  a  constant  stress  plateau  with  values 
close  to  that  at  station  7  in  the  same  region  is  observed. 

The  TJw/(TF  w2)'12  correlation  coefficient  is  positive  at  station  1  and  increases  in  magnitude 
by  station  2  (figure  11).  Station  3  values  above  y+=40  closely  follow  the  station  2  values.  Starting 
with  station  4  the  maximum  positive  value  decreases,  and  near  the  wall  correlation  coefficient  is 
negative.  The  negative  minimum  value  obtained  near  the  wall  decreases  at  stations  downstream 
of  station  4  until  station  6,  and  at  further  downstream  stations  it  increases.  The  location  of  the 

minimum  and  the  location  where  the  uw=0  shift  up  proceeding  downstream.  At  stations  8  and 
9  the  coefficient  changes  its  sign  to  become  negative  in  the  outer  region.  The  positive  peak 

locations  in  the  profiles  downstream  of  station  4  are  further  away  from  the  wall.  By  definition  theuiv 
in  a  2-D  flow  is  zero  due  to  the  flow  symmetry. 

At  the  first  station  uvv/ivt2  is  constant  at  two-thirds  the  level  of  -uv/u?  between  y+=25  and 

400.  At  station  2  this  stress  increases  to  a  level  higher  than  that  of  -Hv/u?.  At  station  3  there  is 

a  slight  decrease  with  respect  to  station  2.  At  station  4  and  other  downstream  stations,  the  values 
change  sign  between  y+=20  to  100.  At  station  6  the  near-wall  maximum  absolute  value  is  5/3 


times  higher  than  that  of  the  maximum  value  of  -uv/u?  at  the  same  station.  At  station  7  the 

absolute  maximum  near  the  wall  decreases.  The  minimum  peak  at  Station  8  is  as  high  as  that 
at  station  7.  At  station  9,  the  near  wall  peak  value  is  decreased  to  =1/2  of  that  of  station  8.  The 

location  where  Uw  changes  sign  shifts  to  higher  y+  values  at  downstream  stations. 

The  vwl(vz'w2y12  correlation  coefficient  is  positive  near  the  wall  at  all  stations  and 
changes  sign  further  away  from  the  wall  proceeding  downstream  (figure  12).  The  minimum  value 
of  the  profiles  decreases  proceeding  downstream  of  station  1  until  station  4  and  is  at  y+=500.  At 
downstream  stations  the  minimum  values  increase;  also  the  near  wall  positive  correlation  region 
extends  farther  away  from  the  wall.  The  maximum  value  is  obtained  at  station  6.  At  downstream 
stations  the  maximum  values  decrease.  The  lateral  pressure  gradient  becomes  negative  at  station 
4  and  stays  negative  at  downstream  stations. 

A  similar  development  is  observed  for  vw/u?.  At  stations  1  through  3  the  value  of  the 

stress  increases.  From  station  4  onwards,  the  near-wall  values  change  sign.  Sign  change  occurs 
at  higher  locations  in  the  downstream  profiles.  The  positive  peak  value  locations  in  the  profiles 
shift  to  higher  y+  values  proceeding  downstream.  The  positive  peak  value  at  station  8  is  =0.47 
which  is  2.5  times  higher  than  at  station  7  and  is  comparable  in  magnitude  to  the  minimum  value 
obtained  at  station  3  which  is  =-0.55. 

The  Uw  and  vw  stresses  have  opposite  signs  to  each  other  and  are  zero  at  locations 

close  to  each  other.  The  major  production  terms  for  Uv,Uw,vw  are  -(UV—+V2—), 

dy  dy 

-(uv^.+vw^li),  respectively.  Production  of  Uw  is  zero  where  -tan  (SSA)=tan 

dy  dy  dy  dy 

(FGA),  or  where  SSA=-FGA.  Production  of  vw  is  zero  where  vw—=-vz^.  Therefore  in  the 

dy  dy 

region  where  SSA=0  and/or  FGA=0,  the  production  of  the  stresses  W,  Uw  and  vw  are 
approximately  zero.  On  the  contrary  -Uv  and  IF  productions  are  always  positive  in  the  layers 
while  IF  production  is  small  and  negative,  except  in  the  outer  layer  where  it  is  small  and  positive. 


CHAPTER  III.  HIGHER  ORDER  PRODUCTS 


The  6  stress  transport  equations  with  boundary  layer  approximations  include  only  the  y 
(vertical)  derivatives  of  the  triple  products  which  are  the  Reynolds  averages  of  the  instantaneous 
Reynolds  stresses  multiplied  with  the  instantaneous  v  fluctuating  velocity,  or  the  flux  of  the 
Reynolds  stresses  in  the  vertical  direction.  Therefore,  the  simplified  equations  do  not  include  the 

1P,W ,  IFw,  uw 2  triple  products.  The  transport  equation  for  the  Reynolds-averaged  triple 
products  can  be  written  as 


DUm.P 

~dT  m 


+PiJk2  +dUk^ijk  eijk 


(4) 


where 

PIJk  1  =m (uJQ  >m+Wk GVQ  >m  ,m 

P«*2=-(W  Uk,m+WK  uj 
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V  •\PlvmUih+Kuku,)j+K  UjUk)j\ 
e;y*=2v(  Ukulm  uJm  +  UjuKm  uim  +  Ujulm  ukm) 

Pfiki  anc*  Pyk2  represent  generation  due  to  Reynolds  stress,  Reynolds  stress  gradient  interaction 
and  mean  strains.  djjk  represents  diffusion,  represents  pressure  interaction,  and  eljk  represents 
dissipation  (Nagano  and  Tagawa,  1 991 ).  ",m"  denotes  differentiation  with  respect  to  that  direction. 
Locations  where  the  Reynolds  stresses  or  triple  products  are  zero  and/or  the  gradients  of  stresses 
or  the  mean  velocity  gradients  are  zero  result  in  reduced  production  for  the  triple  products.  All  the 
triple  products  at  the  wall  are  zero  since  the  fluctuating  velocities  at  the  wall  are  zero.  All  the  triple 
correlations  analyzed  were  found  to  be  of  the  same  order  of  magnitude.  Figures  13  to  22  show 
the  triple  products  nondimensionalized  with  respective  powers  of  the  normal  stresses,  rather  than 
u,3  which  is  not  a  scaling  velocity  for  the  triple  products  as  pointed  out  by  Murlis  et  al.  (1982). 
However  the  discussion  on  the  triple  products  nondimensionalized  with  ut3  is  included  here.  The 
figures  of  the  shear  stresses  and  triple  products  nondimensionalized  with  respective  powers  of 

ut  are  not  presented  in  this  paper.  All  the  triple  products  except  TFw,  v*w,  Ww,  TPshow  similar 


trends  at  both  2-D  and  3-D  stations. 

The  TFv/(lF)  (v2)''2  profiles  show  a  similar  shape  at  every  station  (figure  1 3).  A  minimum 
at  y+=8  and  a  local  maximum  around  y+=40  and  another  local  maximum  around  y+=2000  near  the 
layer  edge  are  observed.  Near  wall  minima  increase  from  1  to  3,  decrease  from  3  to  7,  and 
increase  again  from  7  to  9.  Maxima  near  y+=40  increase  until  station  7  and  decrease  at  further 
downstream  stations.  The  values  decrease  from  their  peak  value  to  form  a  local  minimum  around 
y+=400,  at  stations  5  to  9  these  minima  are  negative  .Compared  to  2-D,  all  3-D  station  values  are 

lower  near  the  wall  up  to  y+=15.  The  y  flux  of  u2  or  TFv  is  contained  in  the  turbulent  kinetic 

energy  (TKE)  equation  and  the  IF  stress  transport  equation  and  (nondimensionalized  with  ut3) 
has  a  local  minimum  at  y+=10  and  a  local  maximum  around  y+=40  at  every  station.  Above  this  y 
location,  the  values  decrease.  This  decrease  is  more  pronounced  progressively  downstream. 
Values  change  sign  between  y+=15  to  20.  Below  this  y+  location  negative  v  fluctuating  velocities 
are  dominant  and  above  it  positive  v  fluctuations  are  dominant.  This  indicates  that  below  y+=15 

to  20,  where  maximum  TKE  production  occurs,  the  I l2  energy  on  the  average  is  carried  towards 

the  wall  and  above  that  y  location  away  from  the  wall.  The  slope  of  IFv/ut3  in  this  region  becomes 

larger  due  to  increased  production  of  IF  with  acceleration  of  the  flow  at  the  last  several  stations. 

The  production  of  IF/2  is  also  the  major  TKE  production  term,  which  means  that  the  TKE 

contained  in  v2 12  and  w2/ 2  mainly  come  from  redistribution  of  IF/2  (Rotta,  1962).  Consequently, 
it  is  expected  that  the  turbulent  diffusion  of  v2  and  w2  will  be  much  smaller  than  that  for  u2. 

The  vw2/  (w2)(v2)1'2  values  are  between  0  to  0.1  for  y+<1000  (figure  14).  This  shows  that 
the  positive  v  fluctuations  on  the  average  occur  with  the  w  fluctuations.  All  3-D  station  values  are 
slightly  lower  than  the  2-D  values.  At  stations  6  and  7  around  y+=20  to  25  negative  minima  are 

observed.  The  y  flux  of  w2  or  vw2  is  small  in  magnitude  with  respect  to  other  triple  products. 
Qualitatively  the  values  increase  towards  the  boundary  layer  edge  and  then  decrease  to  zero. 
Except  at  stations  6  and  7  near  the  wall  all  the  values  are  positive,  indicating  that  w  fluctuations 
occur  more  strongly  together  with  positive  v  fluctuations. 

The  skewness 


shows  the  asymmetry  of  the  probability  distribution  function  (p.d.f)  of  a  fluctuation  quantity.  For 
a  symmetric  distribution  such  as  a  Gaussian  the  skewness  is  zero.  A  positive  skewness  indicates 
a  greater  number  of  occurrences  of  a  quantity  for  large  positive  values  as  compared  to  large 

negative  values.  The  correlation  form  Vs!  ( v2)3'2  which  is  also  the  skewness  of  v  fluctuations 
shows  a  similar  trend  at  all  the  stations  away  from  the  wall  (figure  15).  However,  near  the  wall 
(below  y+=20)  the  variation  is  not  similar.  The  quantity  increases  from  a  negative  value  around 
y+=20  to  a  positive  maximum  around  y+=80  with  higher  values  at  downstream  stations.  Positive 
Sv  indicates  motions  away  from  the  wall  that  are  occasionally  strong.  Around  y+=300  the  values 
are  close  to  each  other  and  are  between  0  and  0.1.  Above  this  height  the  values  monotonically 
increase  till  y+=2000  due  to  outward  motions  and  then  decrease  toward  zero  at  the  boundary  layer 
edge.  For  2-D  flow  near  the  wall  the  values  are  positive  around  =0.2  and  monotonically  decrease 
to  =-0.2  at  y+=20.  The  values  once  more  become  positive  around  y+=40.  The  skewness  Sv  of  v 
fluctuations  below  y+=20  at  stations  1  through  4,  and  at  9  is  negative,  indicating  the  strength  of 
wallward  motions.  At  stations  5  through  8  the  near-wall  Sv  is  slightly  positive.  However,  these 

values  are  relatively  uncertain  because  of  the  small  v*  and  v2  values  in  the  nearest-wall  region. 

The  y  flux  of  v2  or  v2  has  similar  trends  as  vw2  but  with  less  scatter.  Small  negative  values  occur 
very  near  the  wall  due  to  transport  of  v2  toward  the  wall.  In  the  range  y+=40  to  250  a  decrease 
in  the  next  two  stations  is  accompanied  by  an  increase  at  accelerating  downstream  stations  in 
the  same  region.  Values  at  each  station  separately  increase  between  y+=30  and  y+=1000  to  1500; 
near  the  layer  edge  they  decrease.  In  the  outer  layer  the  levels  first  increase  as  the  flow 
decelerates  and  then  decrease  due  to  the  acceleration  at  subsequent  stations. 

The  correlation  uv2/  (IT2)1'2  (v2)  shows  similar  trends  for  all  stations  (figure  1 6).  The  profiles 
form  maxima  at  y+=8  and  decrease  to  values  near  or  below  zero  at  y+=20  to  30,  and  above 
y+=400  they  decrease  till  y+=2000,  and  increase  till  the  layer  edge.  Near  wall  values  are  higher 
than  the  2-D  values.  In  the  y+=30  to  400  range,  while  the  decelerated  flow  stations  have  higher 
values  than  the  2-D,  accelerating  flow  stations  have  lower  values.  The  quantity  shows  that  on  the 
average  the  positive  u  fluctuating  velocities,  whether  they  are  low  in  count  or  not,  result  in  the  flux 

of  the  v2  towards  the  positive  u  fluctuation  direction  near  the  wall.  The  y  flux  of  the  shearing 

stress  uv  or  uv2  (nondimensionalized  with  ut3)  has  a  local  maximum  between  y+=15  to  20  and 
a  progressively  lower  and  negative  local  minimum  between  y+=40  and  60  proceeding  downstream. 
Values  are  close  to  zero  and  positive  and  pass  through  zero  between  y+=200  and  300  at  stations 
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1  through  7,  and  for  the  2-D  profile.  At  stations  8  and  9  values  are  negative  above  y+=20  to  30. 
In  the  outer  layer  the  values  are  negative  at  all  stations  and  the  magnitude  increases,  including 
the  second  station.  Further  downstream  the  magnitude  of  the  negative  values  decreases.  Since 
v2  is  a  positive  definite  quantity,  this  means  that  near  the  wall  below  y+=20  stronger  v  and  stronger 
positive  u  fluctuations  occur  together.  Above  this  y  location  v  fluctuations  are  more  related  to 

negative  u  fluctuations,  resulting  in  negative  uv7  values.  However  above  y+=20  TFv  is  positive 

which  shows  that  u  fluctuations  occur  more  strongly  together  with  positive  v  fluctuations.  These 
observations  taken  together  show  that  in  this  region  stronger  negative  u  fluctuations  and  positive 
stronger  v  fluctuations  occur.  Around  y+=200  the  v  fluctuations  and  ±u  fluctuations  are  evenly 
distributed. 

Negative  v7  and  positive  Uv7  values  near  the  wall  show  a  flux  in  the  average  towards  the 
wall  where  sweeps  are  dominant.  Away  from  the  wall  ejection  motions  play  a  more  important  role. 

At  the  first  3  stations  below  y+=100  the  v7w/(  v7)  (w7)172  values  are  positive,  showing  that 
the  fluctuating  w  velocity  on  the  average  persists  in  the  direction  closer  to  the  tunnel  coordinates 
(figure  17).  At  downstream  stations  especially  at  stations  5,  6  and  7  this  process  results  in 

negative  values  of  V^w/iv7)  (w7)1*.  Above  y+=200  all  station  values  are  close  to  or  below  zero. 
This  is  again  due  to  the  persistence  of  w  fluctuating  velocities  in  the  tunnel  direction.  Near  wall 
positive  values  and  large  negative  values  are  believed  to  be  due  to  the  deflection  of  the  flow  from 
the  direction  of  the  upstream  location  where  it  originated  from  near  the  wall  as  a  result  of  the 
pressure  gradient  being  more  effective  in  this  region.  At  stations  1  through  4  the  y  flux  of  the 

spanwise  shearing  stress  vw  or  vvw  (nondimensionalized  by  uT3)  shows  a  positive  maximum 
below  y+=100  which  is  followed  by  a  negative  maximum  in  the  outer  layers.  At  station  3  the 
positive  maximum  is  at  y+=50  and  negative  maximum  is  at  y+=900.  Proceeding  downstream  until 
station  4  the  near  wall  values  are  negative.  Stations  7,  8,  and  9  show  a  negative  peak  at  y+=30. 

The  v*w  changes  sign  twice  for  station  4  and  at  downstream  stations.  For  stations  1  through  3, 
the  near  wall  v  fluctuations  occur  together  with  stronger  (+)w  fluctuations  and  in  the  outer  layer 
they  occur  simultaneously  with  stronger  (-)w  fluctuations.  For  downstream  stations  where  the  near 
wall  flow  is  turning  back  towards  the  undisturbed  free  stream  flow  direction  the  (-)w  fluctuations 
occur  more  frequently.  At  the  outer  layers  of  different  stations  flow  above  y+=200  to  300  (-)w 
fluctuations  occur  more  frequently  together  with  the  v  fluctuations.  More  (-)w  or  stronger  (-w) 

fluctuations  indicate  that  on  the  average  v7  momentum  is  carried  away  from  the  wing.  However 


at  every  station  there  is  a  y+  range  where  v2w  is  positive. 

The  uvw/u*  is  positive  when  most  of  the  two  of  the  fluctuating  velocities  with  the  same 
order  of  magnitude  have  same  signs  provided  that  the  third  fluctuation  is  positive.  The 

Zmv/ (u7)1'2  (v7)1'2  (w7) 1/2  is  also  positive  when  two  of  the  fluctuating  velocities  are  large  in 
magnitude  and  of  the  same  sign  and  outweigh  the  effect  of  the  opposite  signed  fluctuating 
velocities  in  the  average,  provided  that  the  third  fluctuating  velocity  is  positive.  At  all  stations 

above  y+=400,  uvw/(u7v2w2y,2\s  close  to  or  larger  than  zero  (figure  18).  At  the  first  3  stations 
the  near-wall  values  below  y+=200  are  negative  showing  one  or  all  of  the  fluctuating  velocities  are 
negative  in  the  Reynolds  average  sense.  At  downstream  stations  the  values  are  increasingly 
positive  till  station  6.  At  further  downstream  stations  the  values  decrease  but  they  are  still  positive. 
The  value  of  this  quantity  helps  characterize  the  near-wall  motions.  While  the  w  fluctuating 
velocity  tends  to  stay  in  its  original  direction,  u  and  v  fluctuating  velocities  on  the  average  have 

opposite  signs  to  each  other  (sweep  and  ejection  motions).  At  the  first  3  stations  the  uviv/uT3 
below  y+=200  is  negative,  and  above  this  y  location  it  is  positive  and  increases  downstream  of 
station  1 .  Including  station  4  below  y+=25,  uvw/ uT3  is  positive.  Above  this  y  location  the  values  are 
negative  until  y+=200  to  300.  For  y*  >  300  the  values  are  positive  with  magnitudes  close  to  zero. 

The  IF/ (u2)3'2 shows  the  skewness  (Su)  of  the  u  fluctuating  velocity  (figure  19).  The 
quantity  shows  a  maxima  at  every  station  below  y+=7  due  to  the  wallward  motions,  and  decreases 
towards  zero  at  y  locations  closer  to  the  wall.  The  maximum  value  is  in-between  0.7  and  1 .1 .  The 
values  quickly  decrease  until  y+=40  to  form  a  local  minimum.  A  local  positive  maximum  is  also 
obtained  in  y+=1 00  to  400  region  at  different  y  locations  at  the  various  stations.  Further  away  from 
the  wall  a  local  minimum  at  y+=2000  is  obtained.  Towards  the  layer  edge,  values  increase  towards 
zero.  At  the  first  three  stations  the  minimum  obtained  at  y+=40  is  positive.  At  downstream  stations 
the  values  are  negative.  The  peak  minimum  value  is  obtained  at  station  7  and  decreases  at 
further  downstream  stations.  As  regards  the  2-D  profile,  very  close  to  the  wall,  the  maximum  is 
=0.9  at  y+=4.5,  and  local  minimum  and  local  maximum  obtained  at  y+=40  and  y+=200  are  =-0.1 
and  =0.1  respectively.  High  positive  near  wall  values  are  due  to  wallward  motions  which  carry 
high  speed  fluid  towards  the  wall,  and  high  negative  values  near  the  layer  edge  are  due  to  the 
turbulent  non-turbulent  interface.  The  near  wall  limiting  value  of  Su  for  a  channel  flow  as 
calculated  by  Kim  et  al.  (1987)  using  a  direct  numerical  simulation  (DNS)  is  =0.85.  In  the  current 
study,  very  near  wall  values  show  a  decrease.  This  decrease  was  also  seen  in  the  measurements 


# 


of  Durst,  Jovanovic  and  Sender  (1995).  Gad-el-Hak  and  Bandyopadhyay  (1994)  summarize  the 
Reynolds  number  effects  on  Su  for  2-D  TBLs.  The  discrepancies  found  among  the  measurements 
of  several  researchers  were  attributed  to  the  size  of  the  probe  dimensions.  However, 
measurements  of  Andreopoulos  (1 984)  below  y+  -  6  show  that  Su  values  are  as  high  as  -  1 .7  and 

do  not  change  with  the  Reynolds  number.  The  TF/u*  has  a  positive  peak  at  all  stations  at  y+=10, 
which  denotes  that  most  of  the  u  fluctuating  velocities  are  positive.  Downstream  of  station  3 
where  the  three  dimensionality  of  the  flow  is  more  pronounced,  below  y+=200  the  values  decrease 
to  a  negative  peak  value  at  y+=30.  At  the  third  station  a  local  minimum  is  positive  at  the  same  y+ 

location.  At  y+=200  and  at  stations  1  through  7,  tF/ux3  values  are  closer  to  one  another.  At  the 

last  6  stations  dU/dy  near  the  wall  in  wall  coordinates  is  steeper  than  the  first  three  stations. 
Below  y+=10  very  near  the  wall  the  flow  is  energized  by  the  flow  moving  faster  in  the  region  close 
to  but  away  from  the  wall  which  results  in  positive  u  fluctuations.  Between  y+=10  and  30  the 
energized  flow  very  near  the  wall  lifts  up  and  results  in  higher  negative  fluctuations  until  y+«30  to 
show  a  peak.  The  effect  of  this  lifted  up  low  momentum  fluid  decreases  further  in  the  layers  and 

at  y+«200  this  process  does  not  have  any  effect  in  the  layers.  The  TP  and  TPv  complement  each 

other  in  explaining  this  process.  At  y+=10,  TFv/ux3  reaches  a  minimum  at  all  stations  and  a 
maximum  between  y+=30  to  40  indicating  that  u  fluctuations  are  more  accompanied  with 
increasingly  positive  v  fluctuations  till  this  latter  y  location.  Measurements  of  Murlis,  Tsai  ,and 
Bradshaw  (1982)  in  a  zero  pressure  gradient  2D  TBL  in  a  Ree  range  of  800-4750  in  the  outer 

layers  of  the  flow  show  that  the  TP/u3  triple  product  decreases  with  increasing  Ree.  Their  shear 
stress  and  TKE  balances  in  this  outer  region  showed  that  the  turbulent  transport  mechanism  did 

not  change  qualitatively  with  the  Reynolds  number;  Uv2,  especially,  changed  very  small  in 

streamwise  direction.  They  further  discuss  that  in  this  region  triple  products  are  related  to  shear 
stresses  rather  than  shear  stress  gradients  and  that  turbulent  transport  assumptions  made  in 
calculation  methods  need  not  be  adjusted  for  local  Re9.  However  the  very  large  variations  below 

y+=100  especially  in  TPv  and  TP  suggest  a  closer  look  in  the  near  wall  region. 

The  distribution  of  TPw/(iP)( iv2)1*  shown  in  figure  20  also  supports  the  observations 

made  for  v7w/(v2)(w2)'/2.  At  the  first  four  stations  u^w/ (IF)  (w2)1'2  is  positive  near  the  wall  below 
y+=300.  At  the  downstream  stations  this  quantity  decreases  until  station  7  and  increases  at 
stations  8  and  9.  Above  y+=300,  at  all  stations,  the  values  are  negative  except  for  values  the  near 


18 


layer  edge  at  stations  8  and  9.  The  TFwl ut3  shows  that  at  the  first  4  stations  where  the  three- 
dimensionality  of  the  flow  is  developing,  the  u  fluctuations  occur  accompanied  by  more  (+)w 
fluctuations.  Starting  at  station  5  the  (-)w  fluctuations  are  dominant  below  y+=25.  Above  this  y 
location  mostly  positive  w  fluctuations  accompany  the  u  fluctuations.  Above  y+=400  the  quantity 
is  negative  at  all  stations  with  larger  magnitudes  at  stations  1 ,  2,  and  3  than  the  other  stations. 

Positive  TFw  in  the  y+~30  to  400  range  for  all  the  stations  indicate  that  H7  momentum  is 

associated  with  the  flow  closer  to  the  undisturbed  flow  direction.  For  station  9,  the  values  are 
negative  until  y+=100. 

The  uw'/iw7)  (u7) 1/2 profiles  show  a  maximum  at  y+=7  (figure  21).  The  location  where  the 
profiles  cross  zero  value  changes  from  station  to  station.  At  stations  6  and  downstream  stations, 
the  profiles  cross  zero  at  a  lower  y+  value.  Values  monotonically  decrease  until  y+=2000,  and  then 
increase  toward  the  layer  edge.  This  quantity  shows  that  near  the  wall,  w  fluctuating  velocities 

on  the  average  occur  together  with  stronger  positive  u  fluctuating  velocities.  The  uiv*/ut3  has  a 
positive  peak  at  y+=15  to  20  at  different  stations  and  the  magnitude  of  the  peak  decreases 
proceeding  downstream  until  station  7.  Station  8  and  9  have  peak  values  that  are  close  to  these 
for  station  1 .  Station  1  through  7  values  are  approximately  equal  to  each  other  at  y+=300.  Above 
this  y  location  this  quantity  shows  a  minimum  peak  which  decreases  in  magnitude  proceeding 

downstream  of  station  1  until  7,  and  increases  at  stations  8  and  9.  At  stations  8  and  9UW2 /  u,3 
is  negative  above  y+=40-60. 

The  w*/( w2)312  is  the  skewness  of  w  fluctuations  (figure  22).  Far  upstream,  the  skewness 
is  near  zero,  since  the  mean  flow  is  almost  2D.  For  stations  1  and  2,  the  mean  flow  turns  a  bit 
(table  1)  while  some  turbulence  effects  lag  and  cause  Sw  to  become  positive  in  the  local  wall- 
stress  co-ordinates.  The  skewness  Sw  is  progressively  more  negative  very  near  the  wall  at 
downstream  stations.  The  Sw  in  the  y+  =  100  region  lags  the  inner  region,  i.e.,  Sw  is  still  positive 
until  station  5.  At  station  6  in  the  y+=100  to  1000  range,  the  values  are  close  to  zero.  At  station 
7,  8  and  9  all  the  values  are  below  zero,  except  that  they  are  above  y+=1000  at  stations  8  and 
9.  At  all  stations  in  the  outer  regions,  Sw  decreases  monotonically.  At  the  layer  edge,  the  values 

increase  to  near  zero.  The  TF/ut3  does  not  follow  any  specific  shape.  At  stations  1  through  3 
increasing  positive  values  near  the  wall  are  accompanied  with  decreasing  negative  values  in  the 
outer  layer.  At  station  4  and  at  downstream  stations  the  near  wall  values  are  negative.  At  stations 
4  and  5  in  the  y+=20  to  300  range  the  quantity  is  positive.  Further  decrease  at  station  6  and  at 


downstream  stations  result  in  completely  negative  values  of  the  quantity  which  indicates  that  in 
the  entire  layer  (-)w  fluctuating  velocities  are  dominant. 

Turbulent  diffusion  of  triple  products  includes  derivatives  of  the  fourth  order  products.  The 


triple  products  involved  in  describing  the  TKE  transport  are  7Fv,v*,and,  vw2.  If  the  boundary 

layer  approximations  also  hold  for  the  triple  product  transport  equations,  i.e.,  the  y  variation  is 
much  larger  than  the  x  and  z  variations,  and  the  pressure  interaction  and  dissipation  terms  are 
negligible,  the  following  equations  for  the  triple  products  in  the  TKE  equation  can  be  written: 


TF—  +2uv^  -  (UF}^.+2UF^l)  - 
Dt  By  By  By  By  By 


_ 3-pdV _  Bv* 

Dt  By  By  By 


(6) 


Dt 


By  By  By  By 


Bvzw: 

dy 


Among  the  independent  15  quadruple  products,  the  tFv2,  v*,  and  v ^w2  terms  involved 
in  describing  the  TKE  diffusion  are  presented  as- well  as  the  statistically  important  quantities 
TF and  w*. 

The  !FF/(iF)(lF)  profiles  show  that  in  the  y+=20  to  1500  range  the  distribution  is  flat, 
including  all  the  stations,  with  values  between  1 .0  to  1 .4  with  most  of  the  values  clustered  around 
1.15  (figure  23).  Near-wall  values  are  higher  than  those  for  2-D  with  values  reaching  1 .7  at  station 

7.  Also  near  the  layer  edge,  the  values  increase  monotonically  at  every  station.  The  TFF/u*  at 

the  first  2  stations  profiles  has  approximately  the  same  shape  and  magnitude.  The  profiles  are 
bell  shaped  and  peak  at  y+«  200.  The  peak  values  decrease  until  station  5.  The  decrease  at 
y+=200  is  =  75  %.  At  stations  6  and  7  the  profiles  show  another  peak  at  y+=30  to  40  resulting  in 

double  peaked  profiles.  These  peak  locations  are  close  to  the  locations  where  ~uv/u?,  and 7Fv/u? 

profiles  also  show  peaks.  At  stations  8  and  9  the  values  are  close  to  the  values  obtained  at 
station  1 .  Station  8  and  9  profiles  also  show  double  peaks. 

The  flatness  is  defined  as: 


# 


# 


The  flatness  of  a  Gaussian  p.d.f.  is  3.  If  the  flatness  is  larger  than  3,  the  p.d.f.  is  more  peaked 
than  a  Gaussian  p.d.f.;  if  less  than  3  it  indicates  a  more  flat  top  p.d.f. 

The  V/(F)2  denotes  the  flatness  of  the  v  fluctuating  velocities.  The  profiles  show  that  the 
distribution  of  the  quantity  is  similar  between  y+=50  till  the  layer  edge  (figure  24).  The  values 
decrease  monotonically  from  =3.5  to  =3  until  y+=1000,  and  increase  further  in  the  layers.  Near  the 
wall,  values  as  high  as  =9  are  observed.  The  DNS  calculations  of  Kim  et  al.  (1987)  show  that 
there  is  no  limiting  value  for  this  quantity  at  the  wall.  However,  the  present  data  show  that  very 

near  the  wall  the  values  approach  the  Gaussian  distribution  value.  The  v*/t/,4  values  are  close 
to  zero  near  the  wall  and  at  the  layer  edge.  Profiles  show  monotonically  increasing  values  with 

the  peaks  at  y+=250  to  500.  At  station  2,  "v^/u/  increases  with  respect  to  station  1  values. 
Including  station  5  values  monotonically  decrease.  The  decrease  in  the  peak  value  at  station  5 
is  =70  %  of  station  1  peak  value.  Station  6, 7  values  are  higher  than  those  at  station  5;  with  larger 
values  at  station  6.  Station  8  and  9  values  increase  progressively,  however  the  maximum  value 
lies  in  between  station  3  and  4  values;  which  is  about  60%  of  station  1  maximum  value. 

The  v7wz/(v2)(w2)  values  in  the  range  of  y+=40  to  1 000  decrease  monotonically  and  follow 
each  other  closely  at  every  station  (figure  25).  Near  the  wall  peaks  with  values  as  high  as  =1.65 
are  observed  at  y+=10.  Above  y+=1000,  the  values  increase  monotonically  until  the  layer  edge. 

The  FwV  has  characteristics  similar  to  v*/uf.  Station  2  values  are  higher  than  station  1 

values.  Peak  values  of  the  profiles  decrease  up  to  station  5;  decrease  of  peak  value  at  station 
5  is  =60%  of  the  station  1  value.  At  stations  8  and  9,  the  values  recover  to  80%  of  the  station  1 
value,  and  follow  each  other  closely. 

The  flatness  Fu  of  the  u  fluctuating  velocity  is  large  near  the  wall  and  decreases  to  a  local 
minima  at  y+  =  15  to  20  at  each  station  (figure  26).  At  downstream  stations  it  gradually  increases 
above  y+  =  50  to  near-Gaussian  distribution  values.  The  limiting  value  for  Fu  calculated  by  the 
DNS  method  for  channel  flow  (Kim,  1987)  is  approximately  4.1  at  the  wall.  Near  the  layer  edge, 
above  y+  =  600,  Fu  increases  and  then  decreases  due  to  the  turbulence  intermittency.  The 
flatness  values  for  2-D  TBL  obtained  by  Andreopoulos  (1984)  show  that  Fu  values  very  near  the 
wall  increase  with  decreasing  Ree  value;  values  as  high  as  6.0  were  reported. 

The  limiting  value  for  the  flatness  Fw  of  the  w  fluctuating  velocity  calculated  by  DNS  at  the 
wall  is  approximately  6.75  (Kim,  1987)  (figure  27).  The  data  here  show  near-wall  values  of  4  to 
5,  and  decreases  farther  away  from  the  wall  towards  the  Gaussian  value  of  3.  Near  the  layer 
outer  edge,  the  values  increase  to  about  4  in  the  intermittently  turbulent  region. 
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CHAPTER  IV.  SOME  STRUCTURAL  RELATIONSHIPS 


IV  A.  FLOW  CHARACTERISTICS  IN  THE  VICINITY  OF  W-=0  and/or 

~UW=0 


In  the  vicinity  of  the  locations  where  vw  and/or  uw  change  sign  the  following  observations 
were  made.  The  observations  were  made  in  the  wall-stress  coordinates: 

a)  Derivatives  of  the  quantities  vw  ,7*  ,uv  with  respect  to  ln(y)  are  constants  and 
derivatives  of  uvw  ,  v*w  with  respect  to  ln(y)  are  approximately  zero. 

b)  The  -Vw  and  -Vw  stresses  vary  semi-logarithmically. 

c)  The  stresses  are  not  zero  at  the  same  point  in  the  profiles  but  are  zero  in  the  vicinity 
of  one  another.  The  location  in  the  profiles  where  stresses  are  zero  is  co-ordinate  system 
dependent  (table  4). 

IV  B.  SOME  STRUCTURAL  PARAMETERS  AND  COORDINATE  SYSTEM 
DEPENDENCE 

One  of  the  main  objectives  of  experimental  research  is  to  describe  and  explain  the 
relations  among  flow  variables  in  order  to  reduce  the  independent  variables  in  the  governing 
equations.  In  this  section  some  relations  observed  in  wall  stress  coordinates  will  be  presented. 
The  effect  of  coordinate  system  is  also  discussed. 

Figure  10  shows  the  Uv/(l l2  VZY!Z  stress  distribution  in  wall-stress  co-ordinates.  It  was 

observed  that  the  shear  stress  did  not  scale  on  the  skin  friction  velocity.  Near  the  wall  theTP 
normal  stress  is  observed  as  the  scaling  stress  for  the  longitudinal  shear  stress  component  up 
to  y+=30  (figure  28).  The  Uv/u?  increases  below  y+=30  and  the  correlation  coefficient  also 
increases  towards  zero  value  below  y+=20,  which  indicates  that  the  u  and  v  fluctuating  velocities 
are  not  well  correlated  in  this  region.  However  TF/u?  shows  a  positive  peak  around  y+=10 


indicating  the  presence  of  inactive  motion  in  this  region  (i.e,  high  u  fluctuations  but  low  uv  stress). 

The  quantity  ( — — — )  in  this  co-ordinate  system  denotes  the  ratio  of  the  shear  stress  to  normal 
TF-v2 

stress  difference  (figure  29).  The  ratios  at  each  station  follow  each  other  closely  in  wall-stress 
coordinates  up  to  y+~50.  This  quantity  is  also  the  ratio  of  the  average  momentum  transfer  in  the 
vertical  direction  to  the  difference  in  the  average  stored  kinetic  energy  between  the  horizontal  and 
the  vertical  directions.  This  ratio  shows  that  even  in  wall  stress  coordinates  the  near  wall  structure 

of  Uv  is  governed  by  the  turbulent  interaction  of  fluctuating  velocities. 

Other  quantities  tested  are  the  ratios  of  the  shear  stresses  -Uw/-Uv,  the  ratios  of  the 
shear  stresses  with  respect  to  normal  stresses  -UvFF,  -Uv/w2,  -Uw/u2, 
-Uw/v2,  -Uwfw2,  -vwFF,  -vwFF,  - vwfw 2 ,  and  the  ratio  among  the  normal  stresses,  v2/u2,w2/u2. 

Among  these  quantities,  uv/v2=-0.5  for  all  stations  between  y+*100  to  1000.  The  -vw/v2  is 
constant  in  the  same  range  with  varying  constants  for  stations  1  through  5.  At  stations  6  and  7 
the  stress  distribution  varies  throughout  the  layers.  However  |i/p  \FF  is  shown  to  be  nearly  a 

constant  (Olgmen,  and  Simpson,  1 995a).  All  other  quantities  did  not  show  a  collapse  of  the  data 
on  top  of  another  at  every  station. 

Different  co-ordinate  axes  tested  did  not  result  in  better  presentation  of  the  data.  In  the 
first  co-ordinate  system  used,  the  x  axis  was  aligned  with  the  shear  stress  vector  direction  at  each 
point  of  the  profiles;  the  y  axis  was  kept  perpendicular  to  the  surface  and  z  axis  completed  a 
right-handed  co-ordinate  system  and  all  other  stresses  were  transformed  using  a  tensor 

transformation.  By  definition,  at  each  point  -vw  stress  is  zero.  The  -Uw  stress  and  the  normal 

stresses  did  not  scale  on  the  u2  in  this  coordinate  system  and  Uwlu*  stress  magnitude  was  not 
zero,  but  was  smaller  in  this  co-ordinate  system  compared  to  the  wall  stress  coordinate  values. 
This  was  also  observed  in  simulation  studies  of  Senstad  and  Moin  (1992).  The  -Uvfu2  profiles 
did  not  correlate  as  well  in  this  coordinate  system  as  they  did  in  wall-stress  co-ordinates. 

The  next  co-ordinate  system  was  selected  to  investigate  the  relation  between  the  -vw 
shear  stress  and  the  other  stresses.  For  the  present  flow  studied,  due  to  sign  change  of  the 
lateral  pressure  gradient,  shear  stresses  -Uw  and  -vw,  both  in  free-stream  coordinates  and  wall- 
stress  co-ordinates,  are  observed  to  change  sign  within  the  layers.  It  is  difficult  to  examine 


relations  between  the  normal  stresses  and  the  -vw  shear  stress.  In  the  second  co-ordinate 


system  the  x  axis  was  defined  such  that  the  minimum  value  of  -vw  at  each  station  would  be 
zero.  The  y  axis  was  kept  perpendicular  to  the  surface.  The  tensor  transformation  uses  only  one 
angle  at  each  station.  There  was  not  a  specific  distribution  found  between  the  -vw  stress  and 
any  other  stresses. 

Principal  stresses  and  principal  stress  directions  were  calculated  using  wall  stress  data. 
By  definition  in  principal  stress  coordinates  the  shear  stresses  are  zero.  Principal  stress 
coordinates  are  defined  with  the  unit  vectors  along  the  direction  of  principal  stress  vectors  and 
these  unit  vectors  are  orthogonal  to  each  other.  Principal  stress  coordinates  are  defined 
separately  at  each  y  location  in  the  profiles.  The  maximum  shear  stress  magnitude  in  principal 
stress  coordinates  is  equal  to  the  absolute  value  of  half  of  the  difference  between  two  principal 
stresses.  Its  direction  lies  in  the  plane  defined  by  the  unit  vectors  of  the  principal  stresses  and 
makes  45°  to  each  of  the  unit  vectors.  By  definition  shear  stress  magnitudes  are  positive.  Also 

by  definition  the  diadic  (TFFF+wF)  of  the  Reynolds  stress  tensor  is  a  constant  and  equals  the 


sum  of  the  principal  stresses  (TF~pFF~pFviF~p).  The  profiles  of  the  normal  stresses  in  wall  stress 
coordinates  and  the  principal  stresses  in  principal  stress  coordinates  are  similar  to  each  other. 
The  TF~p/  u?  values  are  =0.5  larger  than  TF/u?  above  y+=60.  The  maximum  values  of ~v7~pl  uf 

profiles  are  =0.3  less  than  IF  I  u?  profiles  and  Wpl  u*  values  are  =0.2  lower  then  W/u?  values. 

No  new  information  could  be  extracted  in  this  coordinate  system. 

Another  structural  parameter  can  be  defined  as 


The  ratio  shows  the  distribution  of  longitudinal  and  spanwise  fluctuations  with  respect  to 
vertical  fluctuations  (figure  30).  In  purely  isotropic  flow  the  ratio  equals  2.  Near  the  wall  S, 
becomes  large  showing  the  anisotropic  nature  of  the  TKE  distribution  in  this  region.  This  is 
essentially  due  to  the  dampening  of  "v"  fluctuations.  In  the  viscous  sublayer  the  ratio  approaches 
infinity  as  the  wall  is  approached  since  u’  and  w’  =  y,  and  v’  =  y2  due  to  continuity  requirements. 
However  the  data  indicate  a  decrease  very  near  the  wall.  At  stations  1  to  3  the  ratio  varies  semi- 
logarithmically  from  y+=100  (  S,=5)  to  1000  (S,=3.5).  At  downstream  stations  where  the  three- 
dimensionality  of  the  flow  is  increasing,  the  ratio  decreases  in  this  y+  range  with  values  higher 


than  2  throughout  the  layers.  The  values  observed  in  this  range  are  comparable  to  the  data  of 
Schwarz  and  Bradshaw  (1992).  In  the  y+=1 00  to  1000  range,  TF/u?  decreases  at  the  downstream 

stations;  however  the  variation  of  v*/uf  and  wF/u?  are  not  monotonical.  Littell  and  Eaton  (1991) 

show  the  parameter  calculated  for  several  3-D  flows  including  Anderson  and  Eaton  (1989), 
Bradshaw  and  Pontikos  (1985),  Driver  and  Johnston  (1990),  Elsenaar  and  Boelsma  (1974),  and 
the  2-D  data  of  Klebanoff  (1955).  Their  figure  4.14  shows  larger  values  of  S1  of  the  order  of  6  to 
10  as  the  wall  is  approached.  The  variation  of  this  ratio  in  the  present  flow  is  closer  to  the 
measurements  of  Anderson  and  Eaton,  and  Klebanoff. 

IV  C.  TRANSPORT  VELOCITY  OF  TURBULENT  KINETIC  ENERGY 

(TKE) 


The  convection  velocity  of  TKE  is  defined  as  an  average  velocity  as  a  means  to  explain 
the  transport  of  the  TKE  as  (Townsend,  1976,  p.125) 

tr+vz+wz 

The  direction  of  the  energy  flux  is  expected  to  be  from  the  regions  of  high  intensity  to  the 
regions  of  low  intensity.  A  similar  definition  was  used  by  Bradshaw  (1967)  in  order  to  model  the 
turbulent  diffusion  term  in  the  transport  equation  of  TKE.  If  such  a  velocity  exists  the  TKE  of  the 
flow  can  be  defined  to  express  the  turbulent  diffusion. 

Continuity  relations  requires  to  be  zero  at  the  wall.  Below  y+  <  20  the  velocity  is 
negative;  the  profiles  show  a  local  minima  between  y+=7  and  10  (figure  31).  Negative  denotes 
that  the  turbulent  diffusion  is  towards  the  wall.  TKE  production  is  maximum  near  the  wall  between 
y+=1 5  to  20,  and  TKE  is  diffused  in  both  directions  around  this  location.  At  the  first  two  stations 
between  y+=40  to  400,  Vq2  is  nearly  a  constant.  At  stations  3  and  4  there  is  a  decrease  around 
y+=1 50.  At  stations  5, 6,  and  7  the  velocity  is  doubled  around  y+=40  indicating  a  larger  penetration 
velocity  of  the  diffusing  TKE.  Profiles  at  these  stations  show  a  minimum  between  y+=200  to  300. 
Towards  the  layer  edge,  at  all  stations,  Vq2  increases;  outside  the  boundary  layer  it  should  go  to 
zero. 

Bradshaw  and  Pontikos’  (1985)  measurements  on  an  "infinite’"'  swept  wing  include  the 
vertical  transport  velocity  above  y+=230.  Data  showed  that  the  velocity  increased  till  the  layer 
edge  with  negative  numbers  at  the  last  2  stations  around  y+=250.  The  “strong  decrease"  of  Vq2 


with  increased  three-dimensionality  observed  for  that  flow  is  not  seen  in  the  present  data. 
Schwarz  and  Bradshaw  (1993)  studied  the  development  and  decay  of  a  3DTBL  generated  by  a 
30°  bend  above  y+=150.  The  data  show  that  in  the  development  stage  Vq2  increases  towards  the 
layer  edge.  With  increased  three-dimensionality  of  the  flow,  is  higher  at  the  layer  edge  and 
gradually  decreases  below  zero  towards  the  wall.  In  the  section  where  the  crossflow  decays,  Vq2 
near  the  wall  gradually  increases  and  becomes  positive.  The  flow  in  both  regions  were  subjected 
to  adverse  pressure  gradient  conditions.  However  for  the  present  measurements  the  streamwise 
pressure  gradient  is  favorable  at  station  2  until  station  6.  The  values  obtained  for  those  flows 
near  the  layer  edge  are  slightly  higher  than  the  values  obtained  for  the  present  flow. 


V.  CONCLUSIONS 


Experimental  results  for  some  higher  order  turbulence  statistics  were  presented  for  a 
three-dimensional  pressure-driven  turbulent  boundary  layer  generated  by  an  idealized  wing-body 
junction.  The  data  included  the  mean  velocity,  Reynolds’  stresses,  10  triple  products  and  selected 
quadruple  products  at  9  stations  and  at  one  mean  two-dimensional  flow  station.  The  triple 
products  indicate  that  the  near  wall  region  is  dominated  by  the  sweep  motions  (+u,  -v),  and  above 
y+=15  to  40  by  the  ejection  motions. 

Turbulent  kinetic  energy  is  transported  towards  the  wall  below  y+=15  to  20  and  is 
transported  away  from  the  wall  above  this  y+  location,  as  shown  by  the  TKE  transport  velocity. 

The  normal  stress  distribution  is  seen  to  be  highly  anisotropic.  The  ratio  (JP+vF)fF  is  an  order 
of  magnitude  different  than  2,  the  value  obtained  for  isotropic  normal  stress  distribution.  The  ratio 
varies  semi-logarithmically  in  the  y+=100  to  1000  range  for  the  first  three  stations  and  for  the  2-D 

station.  The  ratios  Uv/TP  and  2ui ?/(TP-lP)  correlate  the  data  below  y+=50  at  all  stations. 

Triple  products  together  describe  the  flux  of  the  stresses.  Momentum  is  always  positive 
and  is  transferred  with  fluctuating  velocity  components.  If  a  triple  product  is  positive  it  shows  that 

the  normal  stress  is  transferred  with  a  positive  transport  velocity.  The  I P  transfer  is  described  by 
the  TP,  TPv  ,and  TFw  triple  products  in  the  x,  y,  and  z  directions,  respectively.  Large  variations 
of  these  triple  products  occur  mostly  below  y+=200.  The  TP  triple  product  shows  a  positive  peak 
at  y+=10,  while  at  the  same  y+  location  TFv  shows  a  negative  peak,  which  together  indicate  that 
the  TP  momentum  is  carried  towards  the  wall  with  faster  moving  fluid  in  the  xwc  direction.  Above 
y+=10  this  effect  is  decreased.  Between  y+»15  and  20,  TFv  changes  sign  showing  that  thelF 

momentum  is  carried  away  from  the  wall.  The  positive  peak  is  located  at  y+=30  to  40.  ThetP 
triple  product  changes  sign  at  station  4  and  downstream  stations  between  y+=15  and  20,  with  a 
negative  peak  at  y+=30.  At  the  first  4  stations  the  TFw  triple  product  is  positive  till  y+=  200  to  400, 

which  shows  that  the  TP  momentum  transfer  persists  in  its  original  direction.  At  the  downstream 
stations  TFw  is  negative  near  the  wall  indicating  that  the  TP  momentum  is  transported  away  from 


27 


the  wing.  Therefore  very  near  the  wall  below  y+=15  to  20  the  U2  momentum  is  transferred 
towards  the  wall  and  above  this  y+  location  it  is  transferred  away  from  the  wall  with  increased 
transport  velocity.  The  three  dimensionality  of  the  flow  increases  this  transport  velocity.  The  tP 
momentum  transfer  does  not  seem  to  be  affected  by  the  moderate  lateral  pressure  gradients  in 
the  whole  layers.  At  station  5  and  downstream  stations  very  near  the  wall  the  momentum  transfer 
persists  in  its  original  direction  however  with  negative  w  fluctuations  below  y+=30. 

The  v2  flux  is  described  by  the  uvz,  v2  and  vzw  triple  products.  At  all  stations  the  v7 

triple  product  below  y+=  30  is  less  than  zero  indicating  that  the  v2  momentum  is  transferred 

towards  the  wall.  Above  this  y+  location  the  v2  triple  product  is  positive,  indicating  that  the  v2 

momentum  is  in  the  average  carried  away  from  the  wall.  The  uv2  profile  has  a  positive  peak 

around  y+=1 5  including  all  the  stations,  and  the  quantity  changes  sign  between  y+=20  to  40.  The  v7 
momentum  is  transported  with  a  positive  transport  velocity  very  near  the  wall  in  the  xwc  direction, 
and  it  is  transported  with  a  negative  transport  velocity  above  y+=20  to  40.  Above  y+=150  however 

the  uv 2  is  more  negative  in  the  first  three  stations  and  recovers  with  the  increased  three 

dimensionality  of  the  flow.  The  \Fw  is  increasingly  positive  at  the  first  three  stations  up  to  y+=300 

which  shows  a  lag  is  involved  with  the  v2  momentum  transfer  in  the  lateral  direction.  However 

above  this  location  the  vzw  is  negative  and  smaller  proceeding  downstream  of  station  1  which 

shows  that  the  pressure  gradients  are  affecting  the  fluctuating  momentum  more  in  the  outer  layer 
than  in  the  inner  and  logarithmic  layers. 

The  w2  transfer  is  described  by  uw2,  vwz  and  w3  triple  products.  The  w2  momentum 
transfer  in  the  xwc  direction  is  with  a  positive  transport  velocity.  This  transport  velocity  value 
decreases  with  increased  three  dimensionality  of  the  flow.  At  stations  6  and  7  the  value  changes 

sign  at  y+=40.  All  station  uw2  values  are  approximately  equal  to  each  other  at  y+=350,  and  above 
this  y+  location  the  values  are  negative  at  every  station.  The  vwz  triple  product  values  show  that 

the  w 2  momentum  is  also  transferred  away  from  the  wall.  At  stations  6  and  7  the  vwz  values  are 
negative  below  y+=40  indicating  that  the  momentum  transfer  is  towards  the  wall,  above  this  y+ 
location  the  vw 2  is  positive.  The  w2  triple  product  shows  that  at  the  first  station  below  y+=50  and 


at  stations  2  and  3  below  y+=350  the  TF  transfer  in  wall  stress  coordinates  is  with  a  positive 

transport  velocity  in  the  zwc  direction  indicating  that  the  TF  transfer  is  also  directionally  lagging. 
With  the  reversed  lateral  pressure  gradient  between  stations  3  and  4  the  very  near  wall  w 
fluctuating  velocities  at  station  4  become  strongly  negative  however  above  y+=20,  the  TF  is 
positive,  indicating  that  the  mid  layer  flow  is  associated  with  strong  positive  w  fluctuations.  With 
increasing  negative  lateral  pressure  gradients  the  lag  in  the  TF  momentum  transfer  is  clearly 
demonstrated.  At  the  last  three  stations  the  TF  is  negative  throughout  the  layer  except  above 

y+=1000.  At  all  stations  above  y+=350  the  IF  values  are  mostly  less  than  zero,  which  presumably 
indicates  that  the  outer  layer  fluctuating  velocity  structure  responds  faster  to  changes  in  the 
pressure  gradients.  This  is  more  clear  especially  at  stations  1,  2  and  3  that  while  the  TF  is 

increasingly  positive  in  the  inner  layer  the  triple  product  becomes  increasingly  negative  in  the 
outer  layer. 
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TABLE  2.  UNCERTAINTIES  IN  MEASURED  QUANTITIES  WITH  20:1  ODDS 


U/Ux  ±0.075 

TFv  /  (u7)(v2y'2 

±  0.020 

IF/(tF)2 

±0.048 

V/l/t  ±0.026 

IFw  /  (Z77)(iv2)1/2 

±  0.012 

TFv/(u2)3/2(vzy'2 

±0.139 

W IUX  ±0.05 

’vsw/(v2)(wzy/2 

±0.011 

U2w/(LF)3'2(w2y12 

±0.078 

TF/ljf  ±  0.08 

uF/(F)1/2(F) 

±0.015 

FF/(7F)C^) 

±0.039 

F/Ut2  ±  0.029 

uw2  /  (tPy^Cw2) 

±  0.009 

u2vw/(tF)(v^)1/2(vF)1/2 

±0.22 

TF/U2  ±0.037 

vw2/(vsy/2(w2) 

±  0.012 

u2w2/(TF)(w2) 

±0.034 

JUv/(u2V2)1'2  ±0.023 

uvwf(u2y/z(v2y/2 

(W^y®  ±  0.006 

’uv3/(u2y,z(v2)3/2 

±0.073 

-Uw/(u7w2y'2 

±0.021  TF/(uz)3/z 

±  0.024 

iyi71!iv/(tF)1/2('i^)(TF)1/2 

±0.023 

jvwi(vz~w2y12 

±0.018  ^/(V2)3'2 

±  0.052 

-Uvvv?/(iP)1/2(172)1/2("iv2) 

±0.023 

Wliw2)™ 

!  ±  0.022 

uw3/(u7y/2(wz)3J2 

±0.020 

F/(F)2 

±0.084 

- 

IFw/i vz)3'2(w2y12 

±0.052 

^iv*/(F)(TF) 

±0.023 

vw3  /  (v2y/2(wz)3/z 

±0.058 

~vFl(w2)2 

±0.065 
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TABLE  3.  SOME  LENGTH  SCALES  OBTAINED  FROM  LDV  DATA  IN  WALL-STRESS  COORDINATES. 


Stations 


8, 

(mm) 

im 

8, 

(mm) 

8. 

(mm) 

85 

(mm) 

6. 

(mm) 

IBM 

+7.14 

+6.90 

+4.89 

-3.60 

+0.57 

-3.02 

-0.43 

+9.14 

+7.54 

+6.30 

-10.43 

+1.70 

-8.73 

-2.80 

+9.61 

+6.86 

+6.84 

-13.57 

+2.47 

-11.10 

-4.79 

+7.45 

+5.53 

+5.68 

-11.43 

+1.54 

-9.88 

-3.47 

+5.95 

+5.37 

+4.66 

-6.10 

+0.52 

-5.59 

-1.08 

+5.30 

+5.24 

+4.13 

+0.20 

-0.21 

-0.01 

-0.09 

+5.24 

+5.20 

+4.10 

+0.54 

-0.19 

+0.35 

-0.05 

+5.16 

+5.08 

+3.97 

+2.03 

-0.32 

+1.71 

-0.12 

+5.71 

+5.68 

+4.33 

+1.59 

-0.23 

+1.37 

-0.06 

+6.2 

+6.2 

+4.5 

-0.29 

+0.03 

-0.26 

0 

U,=Velocity  magnitude  at  the  layer  edge. 

5=  y  where  \[U*+W*  !Ug= 0.995=  Boundary  Layer  Thickness 


OO 

8,=J(1  ~—)dy=  Streamwise  Displacement  Thickness 

o  Ue 

~r  JTp+W2 

82=|(1  -1 - )dy=  Magnitude  Displacement  Thickness 

LJ 

0  e 

83=J(1  ~—)—dy=  Streamwise  Momentum  Thickness 

0  ^0 

7  w 

84=l(- — )dy=  Lateral  Displacement  Thickness 

0  Ug 

”rw  u 

8S=J — (1- — )dy=  Lateral  Momentum  Thickness 

0  Ug  Llg 
" r  \A/n 

86=J  - — ~dy=  Cross  Product  Momentum  Thickness 

0 

Jw2 

- — dy=  Crossflow  Momentum  Thickness 
U2 
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TABLE  4.  THE  y*  LOCATIONS  WHERE  LATERAL  FLOW  GRADIENT  AND  ~UW  AND 

-VW  STRESSES  ARE  EQUAL  TO  ZERO 

■ 

locations  in  th 

dW/dy=0 

e  profiles  where 

-uw=0 

location 

-vw= 0 

location 

station 

free-stream 

co-ordinates 

normal-stress 

co-ordinates 

wall  stress 
co-ordinates 

wall  stress  co-orc 

Jinates 

1 

32.9 

12 

6.6 

- 

- 

2 

33.4 

9.7 

8.1 

- 

- 

3 

38.2 

28.9 

8.9 

5 

10 

4 

62.1 

26.2 

24.9 

15 

25 

5 

81.4 

56.9 

65.9 

45 

50 

6 

124.0 

97.9 

120.5 

80 

110 

7 

175.6 

153.6 

173.5 

95-100 

140-150 

8 

241.3 

241.1 

247.5 

167.3 

218.5 

9 

340.7 

348 

348.1 

207.5 

321.8 
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Figure  1 .  Wing  shape  and  measurement  locations.  Full  arrows  wall-stress  direction.  Empty  arrows 
tree-stream  direction. 
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Figure  2.  U/uT,  mean  velocity  component  in  wall-stress  coordinates.  The  uncertainties  are  on  the 
order  of  the  symbol  size. 
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Johnston  law-of-the-wall  plot 


Figure  4.  V/uT,  mean  velocity  component  in  wall-stress  coordinates.  The  uncertainties  are  on  the 
order  of  the  symbol  size. 
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Figure  6a.  Flow  angle  calculated  in  wall-stress  coordinates.  The  uncertainties  are  on  the  order 
of  the  symbol  size. 
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Figure  6b.  Flow  gradient  angle  calculated  in  wall-stress  coordinates.  The  uncertainties  are  on  the 
order  of  the  symbol  size. 
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Figure  6c.  Turbulent  shear  stress  angle  calculated  in  wall-stress  coordinates.  The  uncertainties 
are  on  the  order  of  the  symbol  size. 


1 


100  2  345  101  2  345  102  2  345  103  2  3 

y* 

Figure  7.  The  IF/uT2  normal  stress  component  in  wall-stress  coordinates.  The  uncertainties  are 
on  the  order  of  the  symbol  size. 
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Figure  8.  The  ~vllu?  normal  stress  component  in  wall-stress  coordinates.  The  uncertainties  are 
on  the  order  of  the  symbol  size. 
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Figure  9.  The  Tv^/u,2  normal  stress  component  in  wall-stress  coordinates.  The  uncertainties  are 
on  the  order  of  the  symbol  size. 


Figure  1 0.  The  uv/ (u7  y7)1'2  shear  stress  correlation  coefficient  in  wall-stress  coordinates.  Typical 
uncertainty  is  ±0.023. 


Figure  12.  The  vwKv2  w*)1/2  shear  stress  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  1 3.The  TFvl (u2)  (v7)'12  triple  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


Figure  14.  The  'vw7/(w7)  ( v ?)1/2  triple  product  correlation  coefficient  in  wall-stress  coordinates. 
Typical  uncertainty  is  ±0.012. 
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Figure  15.  The  ’v2l{y2f2  triple  product  correlation  coefficient  in  wall-stress  coordinates.  Typical 
uncertainty  is  ±0.052. 
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Figure  1 6.  The  uv2/ (v2)  (u2)'12  triple  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


5 


1 


100  234  101  234  102  234  103  23 

y 

Figure  17.  The  v2 w/ (v2)  (w2)V2  triple  product  correlation  coefficient  in  wall-stress  coordinates. 
Typical  uncertainty  is  ±0.01 1 . 


100  234  101  234  102  234  103  23 

y 


Figure  18.  The  uvwlfjpy12  (v2)'12  (w7)172  triple  product  correlation  coefficient  in  wall-stress 
coordinates.  The  uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  19.  The  t73/(tJ?)3/2  triple  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


Figure  21.  The  Uw2/( w2)  (lF)V2  triple  product  correlation  coefficient  in  wall-stress  coordinates. 
The  uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  22.  The  I F/(w7)3'2  triple  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


48 


•  t 


Figure  23.  The  TFv2 / (u2)  (v7)  fourth  order  product  correlation  coefficient  in  wall-stress 
coordinates.  Typical  uncertainty  is  ±0.039. 


Figure  24.  The  v*/(v*)2  fourth  order  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  25.  The  ’vTw2l(y2)  (w7)  fourth  order  product  correlation  coefficient  in  wall-stress 
coordinates.  The  uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  26.  The  TF/(iF)z  fourth  order  product  correlation  coefficient  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  27.  The  w^/fiv7)2  fourth  order  product  correlation  coefficient  in  wall-stress  coordinates. 
The  uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  28.  The  uv/TF  structural  parameter  calculated  in  wall-wall  stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 
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Figure  29.  The  2uv/(u7-'v2)  structural  parameter  calculated  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


Figure  30.  The  (u2+W)/v2  structural  parameter  calculated  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


Figure  30.  The  (u7+w2)/v2  structural  parameter  calculated  in  wall-stress  coordinates.  The 
uncertainties  are  on  the  order  of  the  symbol  size. 


Figure  31 .  The  Vq2  turbulent  kinetic  energy  transport  velocity.  The  uncertainties  are  on  the  order 
of  the  symbol  size. 


TABLE  1.  UNCERTAINTIES  IN  MEASURED  QUANTITIES  WITH  20:1  ODDS 


u \/ux 

±  0.075 

TFv/u? 

±  0.039 

V/uT 

±  0.026 

IFw/u f 

±  0.066 

W/ut 

±0.05 

v^w/u? 

±0.015 

uV 

±0.08 

UF/u? 

±0.016 

F/u2 

±  0.029 

uw'^/u? 

±  0.059 

W/u? 

±  0.037 

vw'd/u? 

±0.012 

- uvlu f 

±  0.043 

Uvw/u f 

±0.019 

-Uw/u* 

±  0.023 

TF/u? 

±0.189 

-vw/u? 

±0.012 

±0.018 

W/u?  . 

±0.05 
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In  this  Appendix  A  figures  of  some  of  the  flow  parameters  in  wall-stress  coordinates  that  are  not 
included  in  the  body  of  the  report  are  presented.  The  x  coordinate  of  the  wall-stress  coordinates 
is  aligned  in  the  direction  of  the  mean  flow  at  the  first  point  of  the  mean  velocity  profile  and 
changes  from  measurement  station  to  measurement  station.  The  "y"  coordinate  is  perpendicular 
to  the  wall  and  the  "z"  axis  completes  the  right-handed  coordinate  system.  In  the  first  section  of 
the  Appendix  figures  of  mean  velocity,  Reynolds’  stresses  and  triple  products  nondimensionalized 
with  respective  powers  of  Uref  are  presented.  Next  the  transport  velocities  for  the  normal  stresses 
and  the  TKE  in  three  directions  are  presented.  Next,  the  mean  velocity,  Reynolds’  stress  and 
triple  product  deivatives  with  respect  to  the  "y"  and  the  production  terms  for  the  triple  product 
transport  euqations  in  the  vertical  direction  using  these  derivatives  are  presented.  Some  structural 
parameters  are  also  presented  in  this  Appendix. 
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Figures  in  WALL  coordinates.  Transport  velocity  plots  for  all  stations  for  IF,  v2,  w2,  q2  on  top 
of  another.  Where  : 

UU2  =TF/TF,  Vul  =TFv/TF,  Wu2  =TFw/TF, 

UV2 =uv2/~v2,  VV2  = v2/  v7,  WV2  =  v2w/~\F, 

UW2  =TJw2/w2,  VW2  =  vvF/'w 2,  WW2  ^wF/w2, 

Uq2  =TF +TFF +uw2/TF +  v2  +  w2,  Vql  =TFv+  v3 + vwF/TF + v2 + w2 ,  Wq2 = IFw+  Fw+vF/TF + F  +  w2 
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atan(w2/u2)  (u2+w2)/u 


Some  structural  parameters  obtained  using  the  Reynolds’  stresses  and  the  triple  products. 
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In  this  Appendix  the  data  are  presented  in  tunnel  coordinates.  The  “x"  axis  of  the  coordinate 
syst.f 1 T  's  a*'9nec* tunnel  centerline.  The  "y"  axis  is  perpendicular  to  the  tunnel  wall  and 
the  "z”  axis  completes  a  right  handed  coordinate  system.  The  mean  velocity,  Reynolds’  stresses 
and  triple  products  are  nondimensionalized  with  the  relative  powers  of  Uref.  Also  in  this  Appendix 
the  correlation  coefficient  form  of  the  second  and  third  order  velocity  products  are  presented. 


TC 


1 

2 

3 

4 

5 

6 

7 

8 
9 

2D 


Appendix  B:  Tunnel  Coordinates 


100 


V2/uref 


100  234  101  234  102  234  103  23 

y* 


0.0030 


0.0025 


0.0020 


0.001 5 


0.001 0 


0.0005 


0.0000 


100  234  101  234  102  234  103  23 

y* 


Appendix  B:  Tunnel  Coordinates 


hOCOOOvJCTJOi-^OJM 


V2u/Uref 


U3/Uref 


* 


t 


t 


Appendix  B:  Tunnel  Coordinates 


107 


w3/Ur»f 


4.0E— 5 


3.0E-5 


2.0E-5 


1  .OE— 5 


-1  .OE— 5 


-2.0E-5 


-3.0E-5 


-4.0E-5 


100  234  101  234  102  234  103  23 

t 


3.0E-5 
1  .OE— 5 
-1  .OE-5 
-3.0E-5 
-5.0E-5 
-7.0E-5 
-9.0E-5 
-1 .1 E— 4 
-1 .3E-4 


-1 .5E-4 


1 00  234  1 01  234  1 02  234  103  23 

t 


Appendix  B:  Tunnel  Coordinates 


Correlation  coefficient  form  of  the  second  and  third  order  products  in  tunnel  coordinates 
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Appendix  C:  Free- stream  Coordinates 
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In  this  Appendix  the  data  are  presented  in  tree-stream  coordinates  coordinates.  The  "x" 
coordinate  direction  of  the  tree-stream  coordinates  is  in  the  direction  of  the  flow  at  the  outer  edge 
of  the  boundary  layer,  "y"  coordinate  is  perpendicular  to  the  wall  and  the  "z"  axis  completes  a 
right  handed  coordinate  system.  The  mean  velocity,  Reynolds’  stresses  and  triple  products  are 
nondimensionalized  with  the  relative  powers  of  Uref.  Also  in  this  Appendix  the  correlation 
coefficient  form  of  the  second  and  third  order  velocity  products  are  presented. 
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Correlation  coefficient  form  of  the  second  and  third  order  products  in  free-stream  coordinates 
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Appendix  D:  Natural  Coordinates 
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In  this  Appendix  the  data  are  presented  in  natural  coordinates.  The  "x”  coordinate  direction  of  the 
natural  coordinates  is  in  the  direction  of  the  mean  velocity  vector  at  that  point.  The  “x"  coordinate 
direction  varies  at  each  profile  location.  The  "y"  coordinate  is  perpendicular  to  the  "x“  axis  and  lies 
within  the  plane  of  the  "x"  axis  and  the  tunnel  coordinates  "y"  axis.  Therefore  the  "y"  axis  in 
natural  coordinates  also  varies  from  location  to  location.  The  "z"  axis  completes  a  right  handed 
coordinate  system.  The  mean  velocity,  Reynolds’  stresses  and  triple  products  are 
nondimensionalized  with  the  relative  powers  of  Ure(.  Also  in  this  Appendix  the  correlation 
coefficient  form  of  the  second  and  third  order  velocity  products  are  presented. 
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Correlation  coefficient  form  of  the  second  and  third  order  products  in  natural  coordinates. 
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Appendix  E:  Principal 
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In  this  Appendix  the  normal  stresses  are  presented  in  the  principal  stress  corrdinates.  The  normal 
stresses  are  nondimensionalized  with  the  uT.  In  the  principal  stress  coordinates  the  coordinate 
system  is  defined  such  that  all  the  shear  stresses  are  zero. 
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In  principal  stress  coordinates  by  definition  shear  stresses  are  zero.  The  following  shear  stress 
values  are  the  maxiumum  shear  stresses  that  are  calculated  as  half  of  the  sum  of  the  normal 
stresses. 
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